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1.  EXECUTIVE  SUMMARY 


This  report  describes  work  accomplished  under  contract  option  1  of  the  EW  Receiver  and 
Processing  Concepts  Evaluation  Program  (RAPCEval  2).  The  Air  Force  Research  Laboratory 
(AFRL)  awarded  this  contract.  This  branch  of  AFRL  is  located  at  Wright  Patterson  Air  Force 
Base,  in  Da)don,  Ohio.  The  Option  1  task  of  this  program  was  awarded  to  Mercer  Engineering 
Research  Center  (MERC)  on  October  15, 1997.  Work  was  completed  on  the  RAPCEval  1  contract 
in  October,  1996,  and  this  contract  was  initiated  in  November  1996.  This  report  describes  work 
accomplished  on  the  contract  from  December  2,  1997  through  March  31,  1999.  The  project 
includes  five  options  and  has  maximum  duration  of  60  months  from  the  date  of  award.  In  support 
of  the  contract,  a  number  of  activities  and  projects  have  been  initiated.  Some  of  these  are  complete, 
and  many  more  are  continuing.  A  number  of  activities  have  spawned  new  directions  for 
investigative  effort. 

The  Program  Research  Standards  Committee  (PRSComm),  established  at  the  outset  of  the  Basic 
task,  has  continued  to  meet  regularly  and  has  provided  valuable  guidance  and  suggestions  both  as 
to  the  direction  of  students  involved  in  the  effort,  and  the  scope  and  emphasis  of  the  research  efforts 
in  general.  Current  membership  of  this  committee  represents  the  Air  Force  Research  Laboratory, 
Robins  Air  Force  Base,  Mercer  University,  and  Mercer  Engineering  Research  Center.  The  specific 
members  are  listed  m  another  section  of  this  report  (section  4,  Research  Support). 

Three  meetings  of  the  PRSComm  were  held  during  contract  Option  1.  Eight  student  research 
reviews  were  presented  and  were  approved  by  the  PRSComm  in  the  course  of  these  meetings  (see 
section  5,  Project  Activity).  Valuable  discussion  and  suggestions  were  provided  for  directing  and 
focusing  the  students'  work.  In  addition,  suggestions  were  received  for  topics’  reduction  in  scope 
or  better  focus.  Frequently,  unanticipated  or  unknown  resources  and  techniques  were  pointed  out 
for  the  benefit  of  the  students.  One  meeting  involved  a  presentation  by  AFRL  researcher  James 
Stephens,  who  delivered  a  visiting  lecture  with  an  emphasis  on  communication  coimtermeasures, 
with  focus  on  advanced  digital  signal  processing  for  electronic  warfare. 

Students  engaged  in  the  RAPCEval  program  have  continually  enjoyed  firuitful  contact  with 
knowledgeable  personnel  at  the  Air  Force  Research  Laboratory  in  their  respective  areas  of  interest. 
They  have  also  interacted  with  experienced  colleagues  at  Robins  Air  Force  Base,  employees  of 
MERC,  Mercer  University  School  of  Engineering  faculty,  and  various  representatives  of  industry. 

It  should  be  noted  that  software  generated  under  this  contract  is  not  goverranent-owned. 

The  RAPCEval  contract  has  stimulated  gratifying  communication  and  collaborative  research 
effort  among  students,  university  faculty,  MERC  personnel,  personnel  at  the  Air  Force  Research 
Laboratory  and  Warner  Robins  AFB  and  industry.  All  parties  have  expressed  satisfaction  with 
the  contract  results. 
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2.  INTRODUCTION 

These  tasks  specify  requirements  for  analytical  and  research  support  of  in-house  research  at  the 
Air  Force  Research  Laboratory,  AFRL/SNR.  There  is  increasing  sophistication,  quantity,  and 
mobility  of  hostile  radars,  AAM,  SAM,  and  AAA  fire  control  systems.  Electronic  warfare 
receivers  for  radar,  electro-optic,  infi'ared,  ultraviolet  missile  warning  and  electronic 
countermeasures  need  operational  upgrades  to  allow  penetrating  aircraft  acceptable  survivability. 
This  encourages  maintenance  of  in-house  labs  to  support  development,  to  evaluate  concepts,  and 
to  test  new  receivers,  processors,  and  software. 

2.1  EW  Receiver  Effort 

Complex  EW  environments  have  caused  employment  of  numerous  receiving  systems. 
Augmentation  of  in-house  capability  for  evaluation,  novel  concept  development,  and  exploitation 
of  new  technology  is  needed.  Computer-aided  simulation  of  new  systems  and  concepts  can  save 
resources.  New  high  speed  analog-to-digital  (A/D)  converter  technology  may  allow  input 
frequencies  to  be  digitized  in  baseband  before  the  crystal  video  detector,  possibly  allowing  real 
time  digitized  frequency,  pulse  width,  and  pulse  amplitude.  Advancing  materials  technology  for 
infrared  (IR)  /  ultraviolet  (UV)  /  radio  frequency  (RF)  energy  offers  the  possibility  of  augmented 
and  combined  sensors.  Investigation  of  these  materials  is  needed  to  reduce  the  kind  and  number 
of  avionics  needed  in  combat.  Filtering  and  discrimination  advances,  both  in  hardware  and 
software,  may  allow  enhancement  of  fielded  EW  systems. 

2.2  EW  Processing  Effort 

A  modem  EW  system  must  face  an  increasing  number  of  hostile  threats  that  are  multimode. 
Sensors  to  intercept  such  threats  now  include  radar  warning  and  electronic  intelligence  systems. 
Information  from  these  sensors  must  be  processed,  the  threat  identified,  and  appropriate 
countermeasures  initiated  to  coimter  these  threats.  An  augmentation  of  in-house  capability  is 
required  to  evaluate  processor  hardware  and  software,  to  exploit  novel  ideas,  and  to  investigate 
advanced  concepts  such  as  artificial  intelligence  to  determine  the  nature  of  the  threat,  and  what 
countermeasures,  if  any,  to  employ. 

2.3  EW  Exciter  Effort 

Digital  exciters  are  being  developed  to  provide  a  flexible  active  ECM  asset  against  a  wide  variety 
of  modem  threats.  The  need  exists  to  evaluate  the  various  exciter  architectures,  advance  and 
develop  unique  concepts,  and  advance  the  digital  exciter  technology  base.  The  novel  concepts 
and  technologies  must  be  evaluated  for  effectiveness  against  the  proposed  application. 

2.4  EW  Antenna  Effort 

The  role  of  antennas  as  the  "eyes  and  ears"  of  the  sensor  suites  continues  to  make  RF  antenna 
technology  development  vital  to  the  Air  Force  mission.  Airborne  antenna  apertures  of  the  future 
will  be  low  cost,  broadband,  low  radar  cross  section  (RCS),  and  multifunction  in  nature  to  earn 
their  way  onto  platforms  where  space  is  at  a  premium. 
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3.  SCOPE 

The  overall  program  consists  of  a  basic  task  and  four  options  that  are  the  conglomerate  of 

different  work  efforts  and  technologies  within  the  Electronic  Warfare  arena.  Detailed 

descriptions  are  given  as  follows: 

•  Basic  Task  -  The  basic  task  will  provide  the  tasks  necessary  to  analyze  software  and 
hardware  approaches  to  perform  the  exploratory  development  of  EW  technology  in  these 
technology  areas:  radar  hardware,  laser  hardware,  infrared  hardware,  and  ultraviolet 
hardware.  The  task  will  analyze  receiver  and  exciter  technology  to  generate  ECM  signals  to 
improve  ECM  system  performance.  In  addition,  the  scope  of  the  basic  task  will  include 
signal  processing  technology  related  to  the  hardware. 

•  Option  1  -  These  tasks  will  be  those  necessary  to  analyze  receiver  technology  for  application 
to  modem  digital  spectrum  estimation  techniques  in  order  to  improve  EW  /  SIGINT  /  ELINT 
/  IR  /  EO  receiver  performance. 

•  Option  2  -  This  option  consists  of  those  tasks  necessary  to  identify  high  risk  design  areas  for 
an  EW  /  SIGINT  /  ELINT  /  IR  /  EO  hardware  approach,  to  perform  exploratory  design 
assessments  for  selected  functions,  and  to  determine  the  degree  of  parallel  processing 
achievable. 

•  Option  3  -  This  option  is  “reserved.” 

•  Option  4  -  These  tasks  are  those  essential  to  EW  /  SIGINT  /  ELINT  /  IR  /  EO  hardware  and 
signal  processing  including  but  not  limited  to  pulse-deinterleaving,  parametric  extraction,  and 
threat  identification. 
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4.  RESEARCH  SUPPORT 


For  support  of  the  overall  contract,  a  "Program  Research  Standards  Committee"  has  been 
established.  Membership  for  this  committee  was  most  recently  updated  March,  1997.  Current 
members  are: 

•  from  the  Air  Force  Research  Laboratory  in  Dayton, 

•  Mr.  Nicholas  Pequignot  (the  program  manager  for  AFRL) 

•  Mr.  Emil  R.  Martinsek 

•  Mr.  Norman  A.  Toto 

•  Dr.  Duane  A.  Warner 

•  Mr.  Paul  J.  Westcott 

•  from  Warner  Robins  Air  Logistics  Center, 

•  Mr.  Steve  Strawn  (the  program  manager  for  WR-ALC) 

•  Mr.  John  LaVecchia 

•  Mr.  Phil  Oliver 

•  Mr.  Ches  Rehburg 

•  Mr.  Larry  Sheets 

•  from  Mercer  University  and  MERC, 

•  Dr.  Tom  Bass  (the  program  manager  for  MERC) 

•  Dr.  David  Barwick,  (chairman  of  the  standards  committee) 

•  Dr.  Aaron  Collins  (Mercer  U) 

•  Dr.  Behnam  Kamali  (Mercer  U) 

•  Dr.  Paul  MacNeil  (Mercer  U) 

The  EW  Receiver  and  Processing  Concepts  Evaluation  Program  was  awarded  to  the  Mercer 
Engineering  Research  Center  (MERC)  by  WPAFBAVL  under  contract  F09603-93-G-0012-0017. 
This  contract  is  administered  through  WR-ALC.  The  overall  program  has  a  funding  ceiling  of 
$499,940.  Incremental  funding  will  be  accomplished  via  a  series  of  contract  options.  The  basic 
contract  is  $99,998.00,  option  (1)  is  $99,998.00,  option  (2)  is  $99,998.00,  option  (3)  is 
“reserved”,  and  option  (4)  is  $50,000.00. 

Contracts  have  been  awarded  for  the  Basic  program.  Option  1,  Option  2,  and  Option  4. 
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5.  PROJECT  ACTIVITY 


5.1  Steering  Committee,  April  1998 

5.1.1  Meeting  Minutes 


EW  Receiver  and  Processing  Concepts  Evaluation  Program 
Program  Research  Standards  Committee  Meeting 
_ Minutes  2  April  1998 _ 


A  meeting  of  the  Program  Research  Standards  Committee  (PRSC)  for  the  EW  Receiver 
and  Processing  Concepts  Evaluation  Program  (RAPCEval2)  was  hosted  by  Mercer  Engineering 
Research  Center  at  8:30  a.m.  in  the  new  MERC  auditorium.  Committee  members  present  were 
Dave  Barwick,  Tom  Bass,  Aaron  Collins,  Benham  Kamali,  Nick  Pequignot,  Norman  Toto, 
Duane  Warner,  Paul  MacNeil,  Phil  Oliver,  and  Ches  Rehberg.  Also  present  were  several 
representatives  from  Mercer  University,  students  scheduled  to  speak,  and  other  students 
considering  participation  in  the  RAPCEval  program.  Five  student  presentations  were  given. 

Four  were  reports  on  research  in  progress,  while  one  was  a  research  proposal. 

After  welcoming  remarks  by  Dr.  Barwick,  Dr.  Tom  Bass  gave  a  brief  description  of  the 
RAPCEval  program  and  introduced  the  students  who  were  scheduled  to  speak.  The  students  and 
their  topics  are  listed  in  the  table. 


Ron  Brinkley 

“Burst  Error  Correction  with  Reed- 
Solomon  Codes” 

Dennis  Ludwig 

“Identification  of  IR  Images  Using  Neural 
Networks” 

Henderson  Benjamin 

“Selection  of  Reed  Solomon  Codes  Using 
Neural  Networks” 

Randy  Ford 

“Passive  Location  via  Evolutionary 

Genetic  Algorithms” 

Dennis  Ludwig 

“Identification  of  IR  Images  Using  Neural 
Networks” 

ESSSSESSSHH! 

“Hardware  Exploitation  of  RAD  Filtration” 

Ron  Brinkley  showed  overheads  introducing  the  concept  of  Reed-Soloman  error 
correction  methods,  pointing  out  that  RS  codes  are  block  codes.  He  also  gave  a  summary  of  the 
types  of  communication  errors,  together  with  interleaving  techniques  designed  to  combat  the 
errors.  Ron  discussed  his  analysis  and  tabulation  of  the  Reed-Soloman  error-correcting  codes  for 
specific  classes  and  burst  error  categories.  Announcement  was  made  of  Ron’s  presentation  on 
this  subject  to  the  IEEE  Pacific  Rim  Conference  (PACRIM  ’97). 


5 


Henderson  Benjamin  reviewed  his  topic  on  the  application  of  neural  networks  to  the 
selection  of  Reed-Solomon  (RS)  codes.  Henderson  showed  many  of  the  details  of  the  manner  in 
which  RS  codes  are  selected  by  a  neural  net  (NN).  Included  were  training  techniques  for  the  NN 
and  performance  statistics  for  the  trained  network.  Henderson  acknowledged  dependence  on  the 
work  of  Ron  Brinkley  for  availability  of  code  selection  tables  that  were  required  for  the  NN 
application.  This  work  is  nearly  complete  and  the  thesis  is  under  final  review.  An  abbreviated 
version  of  the  thesis  is  available  in  the  RAPCEval  annual  report,  dated  3  February  1998. 

Dennis  Ludwig  reviewed  the  pose  estimation  problem  for  an  object  identified  by  an 
infrared  (or  Other)  image.  He  had  planned  to  extend  a  two-dimensional  solution  developed  at  the 
Wright  Laboratory  to  a  general  solution  in  three  dimensions.  The  chief  difficulty  encountered  is 
the  necessity  of  t^ng  available  two-dimensional  images,  and  somehow  recreating  the  third 
dimension.  A  number  of  questions  were  evoked  by  this  presentation,  ranging  from  whether 
three-dimensional  data  is  available,  to  the  specifics  of  how  the  images  are  brought  into 
alignment. 

Tracy  Tillman  presented  additional  details  on  the  development  of  a  imique  method  of 
table  generation  and  table  lookup  for  use  in  RAD  parameter  production.  He  has  discovered  some 
limitations  for  the  technique,  and  plans  to  finish  a  report  on  Ws  discoveries,  calculations,  and 
algorithmic  development  in  the  next  quarter  or  two. 

Dr.  Barwick  chaired  the  committee  discussions  that  followed  the  last  presentation. 
Suggestions  were  gathered  regarding  the  student  presentations.  The  group  agreed  that  the 
communication  research  projects  directed  by  Dr.  Kamali  were  very  pertinent  to  current  AF  needs 
and  should  be  continued.  Several  sensed  that  Dennis  Ludwig  seemed  to  have  lost  focus  of  the 
original  goals,  probably  due  to  the  departure  of  his  thesis  advisor  recently.  Dr.  Bass,  in 
conjunction  with  Dr.  Collins,  will  be  giving  increased  direction  and  steer  him  back  on  track.  A 
question  arose  as  to  whether  the  RAD  lookup  tables  were  in  use  in  the  PLAID  program. 
Currently,  they  are  not,  but  further  evaluation  is  planned  during  the  next  phase  of  R&D.  The 
comment  was  made  that  Randy  Ford’s  schedule  seemed  optimistic.  Another  question  concerned 
the  potential  to  test  the  differential  evolution  solution  technique  on  actual  PLAID  data.  This  will 
be  investigated. 

During  final  discussions.  Norm  Toto  said  that  topics  formerly  being  investigated  at  AFIT 
were  “going  away.”  He  thought  MERC  might  be  able  to  continue  some  of  these  via  RAPCEval. 

Norm  also  suggested  potential  RAPCEval  research  investigations  into  GPS  anti-jamming 
techniques  and  related  technologies. 

At  this  point,  the  meeting  was  adjourned. 
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5.1.2  Meeting  Agenda 


Agenda 

RAPCEval 

program 


STEERING  COMMITTEE 
MEETING 

98ApriI  2  ---  8:30  AMto2:30  PM 
Main  Auditorium,  Mercer  Engineering  Research  Center 
135  Osigian  Boulevard,  Warner  Robins,  GA 


~i 

Meeting 

Nicholas  Pequignot,  AFRL/SNRP,  Facilitator:  Dr.  Tom  Bass 

called  by: 

Air  Force  Research  Laboratory  Program 

Manager 

Committee 

AF  Rsch  Laboratory 

Robins  AFB 

Mercer  University 

Members: 

Mr  Nicholas  Pequignot 

Mr  Steve  Strawn 

Dr  David  Barwick 

Mr  Emil  R  Martinsek 

Mr  John  LaVecchia 

Dr  Tom  Bass 

Mr  Norman  A  Toto 

Mr  Phil  Oliver 

Dr  Aaron  Collins 

Dr  Duane  A  Warner 

Mr  Ches  Rehberg 

Dr  Benham  Kamali 

Mr  Paul  J  Westcott 

Mr.  Larry  Sheets 

Dr  Paul  MacNeil 

Email 

pequigna@sensors.wpafb.af.mil 

inerr@ewir-wr.robins.afmil 

dbarwick@merc.mercer.edu 

martinser@sensors.wpafb. af.mil 

lavecchia@ewir-wr.robins.af.mil 

tbass@merc.mercer.edu 

addresses: 

totona@sensors.wpafb-af.mil 

oliver@ewir-wr.robins.af.mil 

collins_as@mercer.edu 

wamerda@sensors.wpafb.af.miI 

rehberg@ewir-wr.robins.af.mil 

macneil_pe@mercer.peachnet.edu 

westcopj@sensors.wpafb.af.mil 

sheets@ewir-wr.robins.afmiI 

kamali_b@mercer.peachnet.edu 
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Meeting  Agenda,  continued... 


Schedule 

Continental  Breakfast  /  Greetings 

Dr  David  Barwick 

8:30  AM -9:00  AM 

Meeting  Overview 

Dr  Tom  Bass 

9:00  AM  -  9:20  AM 

Student  Proposals 

Dr  Tom  Bass  /  Students 

Burst  Error  Correction  with 
Reed-SPlornOn  Codes 

ROn  Brinkley 

9:20  AM  -  9:55  AM 

Selection  of  Reed  Solomon 

Codes  Using  Neural  Networks 

Henderson  Benjamin 

9:55  AM  -  10:30  AM 

Break 

10:30  AM -10:50  AM 

Student  Proposals  (continued) 

Dr  Tom  Bass  /  Students 

Passive  Location  via  Evolution- 
arv  Crerietic  Alfforithnis 

Randy  Ford 

10:50  AM  - 11:25  AM 

dJ-  y  VJwXXVtlV  i^Jl^V/J.XULJLlXXO. 

Identification  Of  IR  Images  Us¬ 

Mr  Dennis  Ludvvdg 

1 1:25  AM  -  NOON 

ing  Neural  Networks 

Lunch 

NOON -1:00  PM 

Student  Proposals  (continued) 

Dr  Tom  Bass  /  Students 

Hardware  Exploitation  of  RAD 

Mr.  Tracy  Tillman 

1:00  PM-  1:35  PM 

Filtration 

Discussions  &  New  Business 

Dr.  Tom  Bass 

1:35  PM -2:30  PM 

Adjourn 

2:30  PM 
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5.1.3  Attendance  Roster 

The  attendees  at  this  meeting  are  listed  here: 
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Tehan,  Jack  MERC  (North)  (803)  795-7030  jtehan@aol.com 

Toto,  Norman  A.  AFRL/SNRW  (937)  255-4933  toto@sensors.wpafb.af.mil 

Warner,  Duane  AFRL/SUJO  (937)255-4174  warnerda@aa.wpafb.af.mil 


5.1 .4  Overview  of  the  Program  (Dr.  Bass) 

The  Overview  Briefing  of  the  RAPCEval  Program  as  presented  at  this  meeting  is  reproduced  on 
the  next  7  pages. 
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April  2, 1998 


several  (5)  research  projects  are  ongoing 

★  All  research  has  been  determined  by  the  committe 
to  be  useful  to  the  Air  Force,  as  well  as  having  me 
for  graduate  research 
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Tracy  Tillman,  “Hardware  Exploitation  of  RAD 
Filtration” 
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Improved  Communication  of  Time  and  Time 
Measurements”,  unclassified  report  available  from  MERC 
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★  Kirk  Wright,  “Object  Oriented  Modeling  of  the 
AN/ALQ-172”,  publication  in  process 
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Networks 

Tracy  Tillman,  “Hardware  Exploitation  of  RAD  Filtration 


5.1.5  Presentation  By  Ron  Brinkley 

The  student  briefing  presented  by  Ron  Brinkley  at  this  meeting  is  reproduced  on  the  next  32 
pages. 
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Research  Topic:  BURST  ERROR  CORRECTION  WITH  REED-SOLOMON  CODES 
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various  Reed-Solomon  codes.  Also,  the  parameter  titled  burst  error  correction  efficiency  will 
be  introduced.  Burst  error  correction  efficiency  is  a  measured  value  that  will  allow  a  code 
designer  additional  information  for  evaluating  possible  Reed-Solomon  codes.  The  concepts 
of  burst  error  correcting  capability  and  burst  error  correction  efficiency  will  be  expanded  to 
include  block  and  convolutional  interleaving  techniques. 
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[Ski]  B.  Sklar,  Digital  Communication,  Fundamentals  and  Applications,  Englewood  Cliffs,  NJ:  Prentice  Hall,  1988. 
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Cluster  of  Individual  Errors 
Channel  Memory 
Long  Single  Burst 
Multiple  Short  Burst 
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BLOCK  INTERLEAVER 


Interleaver 


Input  Stream:  Al,  A2,  A3,  A4,  A5,  A6,  A7,  A8,  A9,  AlO,  AH,  A12,  Bl,  B2,  B3,  B4,  B5,  B6, 


.,F9,F10,F11,F12 


◄ -  n  symbols  — — - —  . ► 


A1 

A2 

A3 

A4 

A5 

A6 

A7 

A8 

A9 

AlO 

All 

A12 

B1 

B2 

B3 

B4 

B5 

B6 

B7 

B8 

B9 

BIO 

Bll 

B12 

Cl 

C2 

C3 

C4 

C5 

C6 

C7 

C8 

C9 

CIO 

Cll 

C12 

D1 

D3 

D5 

D6 

D7 

D8 

D9 

DIO 

Dll 

D12 

El 

E3 

~ir 

E5 

E6 

E7 

E8 

E9 

ElO 

Ell 

E12 

FI 

"fT 

F3 

F5 

F6 

F6 

F8 

F9 

FIO 

Fll 

F12 

Output  Stream:  E2,  F2  A3,  B3,  C3,  D3,  E3,  F3,  A4, . ,  C 12,  D12,  E12,  F12 


burst 


De-interleaver 


Input  Stream: 


Al3i>  Cl,  Dl^  El,  FI '  A2,  fi2,  C2,  P2,  E2,  F2, 


A3,  B3,  C3,  D3,  E3,  F3,  A4, 


4 -  n  symbols  - ► 


Al 

A2 

A3 

A4 

A5 

A6 

A7 

A8 

A9 

AlO 

All 

A12 

El 

B2 

B3 

B4 

B5 

B6 

B7 

B8 

B9 

BIO 

Bll 

B12 

eii 

;c2 

C3 

C4 

C5 

C6 

C7 

C8 

C9 

CIO 

Cll 

C12 

j>i 

D2 

D3 

D4 

D5 

D6 

D7 

D8 

D9 

DIO 

Dll 

D12 

El 

m 

E3 

E4 

E5 

E6 

E7 

E8 

E9 

ElO 

Ell 

E12 

F3 

F4 

F5 

F6 

F6 

F8 

F9 

FIO 

Fll 

F12 

,C12,D12,E12,F12 


Output  Stream: 


B3,  B4,  B5,  B6, . ,F9,  FIO,  Fll,  F12 
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RF  RECEIVER  AND  PROCESSING 
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[Fori]  G.  D.  Fomey,  Jr.,  “Burst  Correcting  Codes  for  the  Classic  Bursty  Channel,” 

IEEE  Transactions  on  Communication,  Volume  COM-19,  pp.  772-781,  October  1971. 
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[Kam]  B.  Kamali,  Lecture  notes  on  Reed-Solomon  Codes,  EEE  691.80  Spring  Qrt.  1995,  Mercer  University. 
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EXAMPLE  PROBLEM 
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[Cap]  V.  Cappellini,  Data  compression  and  Error  Control  Techniques  with 
Applications,  Academic  Press,  1985. 

[Gall]  R.  G.  Gallager,  Information  Theory  and  Reliable  Communication,  New  York: 
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5.1.6  Presentation  by  Henderson  Benjamin 

The  student  briefing  presented  by  Henderson  Benjamin  at  this  meeting  is  reproduced  on  the  next 
24  pages. 
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Research  Topic;  NEURAL  NETWORK  SYSTEM  THA  T  SELECTS  REED-SOLOMON 
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Collect/Analyze  data  for  shortened  RS  codes  of  length  255 
Collect/Analyze  additional  data  deemed  useful  to  a  designer 
Develop  Neural  Network 
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Jamming,  and  recording  systems. 
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in  interleaved  cases 
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The  trade  off  of  shortened  codes  isfhe  efficiency  of  the  coderate  10  and  12% 
respectively  in  the  above  case. 
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slower  than  semiconductor  devices. 
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^4.40%  Correct  for  the  code  dimension 


5.1.7  Presentation  by  Randy  Ford 


The  student  briefing  presented  by  Randy  Ford  at  this  meeting  is  reproduced  on  the  next  26  pages 
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Research  Topic: 

COMPARISION  OF  DIFFERENTIAL  EVOLUTION  TO  THE  SIMPLEX  METHOD  IN 
OPTIMIZATION  DURING  PASSIVE  EMITTER  LOCATION 
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minimization  technique  for  real-valued  problems  which  is  based  on  genetic 
algorithms.  This  project  will  compare  the  performance  of  this  method,  with  the 
performance  of  the  Nelder-Mead  simplex  method.  Both  methods  will  be 
implemented  in  C++  and  called  from  an  existing  MATLAB  implementation. 
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Which  one  produces  fastest  solution? 


^  S  oo 
O  W)  o^ 
[i,  O  Cn 

S  “  a 
g  ia  < 

«  w  ol 


'-'  o 

W  o  ^ 

^  5  « 

S  z  U 
^  2  ^ 
^  5  ^ 

U  ^  w 

M 

z  hJ 

o 

Pi  > 

H  M 


A'5« 

\%v 


[£)  ^ 

K  53 
53 

li 

|g 

S  ^  . 

^11 

^  ^  ^ 

^  ^  i 

SIS 

f>s  5!  «< 

gfeg 

Is  fe 

Q  ^ 

t,  g  bi 

®  S 

S 

&!  ^ 

P§ 

U  ^ 


^  + 
Qh  + 

e  ^ 


•'^  d  ^ 

>1  a  ^ 

^  2  S 

I— I  ‘ip  I 

<  .H  ^ 

w)  S  ^ 
a  *S  't:^ 


c5 

N  22 
C2 

B  ^ 

.  td 

c  ^ 

4-> 

s  g 

w  g 

Q  § 


-1— ‘ 

£ 

p 

P 

4-» 

(D 

a 

<D 

a 

e 

(U 

1— ( 

a 

'E, 

a 

Ph 

h— i 

HH 

s 

• 

• 

81 


Identical  Platfornis  -  Pentium  PC 

Determine  best  method  based  on  speed  of 
detection. 
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Consideration  will  also  be  given  to  the  frequency  with 
which  a  method  produces  an  unacceptably  slow  solution 
or  no  solution. 
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.  Begin  coding  C++  Nelder-Mead  Implementation 
Semester  2  (Summer  98-1  hour) 

•  Finish  coding  C++  Nelder-Mead  Implementation 
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Finish  coding  C++  DE  Implementation 
Test  and  Compare  Implementations 
Final  Project  Report 
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•  Get  x,y  ellipse  by  projecting  the  CR  ellipsoid  onto  the 
x,y  plane. 
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★  Passive  Emitter  Location  discussion  is  based  on  An 
Efficient  Method  of  Passive  Emitter  Location”  by  Klaus 
Becker  in  the  October  1 992  IEEE  Transactions  on 
Aerospace  and  Electronic  Systems. 
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line  segments  and  polygonal  faces. 

Initialization:  Start  with  N  +  1  initial  points  which  give 
an  initial  simplex. 
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Reflect  and  expand  the  simplex  away  from  the  high  point 
Contract  the  simplex  in  one  dimension  from  the  high  point 
Contract  the  simplex  along  all  dimensions  toward  a  low  point. 
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■  Stop  when  the  decrease  in  function  value  during  a  step  is  less 
than  some  tolerance 

Restart  at  minimum  . 
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William  H.  Press,  Saul  A.  Teukolsky,  William  T. 
Vetterling  and  Brian  P.  Flannery. 
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Algorithm  that  is  especially  useful  m  minimizing 
functions  with  real  variables  and  many  local  minima. 

Differential  Evolution  is  an  Evolutionary  Strategy 
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them  into  a  special  alphabet. 

Usually  use  a  deterministic  selection  scheme  during  crossover. 

Use  additiop  instead  of  symbol-swapping  during  mutation. 
This  makes  incremental  search  easier. 
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member.  Store  these  values  in  a  cost  array. 

’'s. 

■  Each  member  of  the  current  population  takes  its  turn  as  the 
target  vector  for  mutation,  crossover,  and  selection 
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Determine  the  difference  between  two  randomly 
population  members,  x^  - 
Multiply  by  the  scaling  factor  F.  F*(  x^  -  x^,) 
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■  Produce  the  trial  vector  Xj  by  combining  the  target  with  x^.’. 
.  Conduct  D-1  experiments  starting  at  a  randomly  selected 
parameter. 
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generation  array.  If  the  trial  vector  wins,  its  cost  replaces  its 
parent’s  cost  in  the  cost  array. 
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ermination: 

>  Trials  continue  for  a  pre-set  number  of  iterations  or  until  the 
cost  of  the  best  vector  is  less  than  or  equal  to  a  pre-set  value. 
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Next  Generation 
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Test  and  compare  implementations 


5.1.8  Presentation  by  Dennis  Ludwig 


The  student  briefing  presented  by  Dennis  Ludwig  at  this  meeting  is  reproduced  on  the  next  1 2 
pages. 
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RF  RECEIVER  AND  PROCESSING 
CONCEITS  EVALUATION  PROGRAM 
( RAPCEval ) 


Dennis  A.  Ludwig 
MSEE  Program 
UPDATE 


MPRO-l 


RESEARCH  PROPOSAL  PRESENTATION 
Dennis  A.  Ludwig 

WR-ALC/LNERB 
ALQ-155  System  Engineer 


Background  and  Experience: 

Education:  BS  EE  from  Louisiana  Tech,  1981 

MSA  from  Georgia  College,  1988 
Pursuing  MSEE  from  Mercer,  45  hrs.  completed 
WR/ALC:  Electronic/Software  Engineer 


Research  Topic: 

Pose  estimation  using  a  Neural  Network 


RF  RECEIVER  AND  PROCESSING 
CONCEPTS  EVALUATION  PROGRAM 
( RAPCEval ) 


Dennis  A.  Ludwig 
MSEE  Program 
UPDATE 


M  PRO-3 


PROBLEM  STATEMENT 

-*■  Design  a  Neural  Network  system  to  resolve  the  Pose 
Estimation  Problem  of  a  target/  object  in  three 
dimensions.  That  is  the  estimation  of  location  and 
orientation  of  a  detected  target  in  a  Pose  volume. 

*  Given  an  image  of  an  object  in  a  background,  and  a 
model  that  we  would  like  to  match  to  the  object, 
what  is  the  best  pose  for  the  model?  (Ghee  and 
Doria) 
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RF  RECEIVER  AND  PROCESSING 
CONCEPTS  EVALUATION  PROGRAM 
( RAPCEval ) 


PROPOSED  RESEARCH  OBJECTIVE 

Existing:  The  pose  estimation  problem  has  been  solved  for  a 
two-dimensional  image,  involving  three  degrees  of 
freedom  (two  associated  widi  position,  one  with 
orientation). 

Improvement:  To  add  a  third  dimension,  which  adds  three 
parameters  (Z,  Phi,  Psi)  to  the  current  two-dimensional 
Pose  Estimation  NeuraJ  Network  System.  This  will 
correspond  to  six  degrees  of  freedom. 


Dennis  A.  Ludwig 
MSEE  Program 
UPDATE 
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RF  RECEIVER  AND  PROCESSING 
CONCEPTS  EVALUATION  PROGRAM 
( RAPCEval ) 


Dennis  A.  Ludwig 
MSEE  Program 
UPDATE 


TECHNICAL  SECTION 

3-D  Image  Representation 

*  Three  degrees  of  freedom  for  translation  in  the 

X,  Y  and  Z  directions 

«  Three  degrees  of  freedom  for  rotation  or  orientation 
about  the  or  Z  axis  {<£>,'¥  ,  ©  )  respectively. 

•  Scale  will  be  incorporated  into  a  distance. 

★  Note  that  3-D  Pose  volume  has  6  parameters  describing  it. 


TECHNICAL  SECTION 


RF  RECEIVER  AND  PROCESSING 
CONCEPTS  EVALUATION  PROGRAM 
(RAPCEval) 


Dennis  A.  Ludwig 
MSEE  Program 
UPDATE 


System  Overview 


Image 

(Data) 


Minimum  Distance 
Vector  Calculation 


Pose  Estimation 
Neural  Network 


Model  1 
- z - ' 

/ 

'  I 

1 - — 

iPose  Adjustment  | 

1  y 
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RF  RECEIVER  AND  PROCESSING 
CONCEPTS  EVALUATION  PROGRAM 
( RAPCEval ) 


TECHNICAL  SECTION 


Dennis  A.  Ludwig 
MSEE  Program 
UPDATE 


MPRO-9 


Model 

•  An  imitation  of  an  existing  object 

•Contains  a  priori  knowledge  of  the  feature 

•Is  encoded  in  traditional  data  structure  (two  or 
three  dimensional  array) 

•Referred  to  as  a  synthetic  image 


Data  structure  is  two 
dimensional  array  of 


TECHNICAL  SECTION 

Neural  Network 


★  Architecture  (MLFN)  (MLP)  (Backpropagation) 

★  Type  of  Learning  Procedure  (Supervised) 

★  A  set  of  training  patterns  is  used  to  teach  a  NN. 

.  Input/ Output  vectors 
•  Adjust  weights  based  on  output  error 

★  Each  parameter  variation  will  have  an  output  node 
corresponding  to  it. 

★  (AX'AY/AZ'AO/A'J^'A©) 
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TECHNICAL  SECTION:  NEURAL  NETWORK  CONT. 

A  typical  weight  file  for  Five  input  nodes,  Two  hidden  nodes,  and  Six  outpu^ 
nodes; 

0(1,0)  0(1,1)  0(1,2)  0(1,3) 

0(2,0)  0(2,1)  0(2,2)  0(23) 

0(3,0)  0(3,1)  0(3,2)  0(33) 

0(4,0)  0(4,1)  0(4,2)  0(4,3) 

0(5,0)  0(5,1)  0(5,2)  0(53) 

0(6,0)  0(6,1)  0(6,2)  0(63) 

H(0,0)  H(0,1)  H(0,2)  H(0,3)  H(0,4)  H(0,5) 

H(1,0)  H(l,l)  H(l,2)  H(l,3)  H(l,4)  H(l,5) 
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RF  RECEIVER  AND  PROCESSING 
CONCEPTS  EVALUATION  PROGRAM 
(RAPCEval) 


Dennis  A.  Ludwig 
MSEE  Program 
UPDATE 


TECHNICAL  SECTION;  NEURAL  NETWORK  CONT. 


Parameters: 

★  Learning  Rate,  can  be  different  for  each  layer 

★  Momentum,  can  be  different  value  for  each  layer 

★  Activation  Function  Temperature 

★  Bias 

Number  of  patterns  in  Epoch. 


RF  RECEIVER  AND  PROCESSING 
CONCEPTS  EVALUATION  PROGRAM 
( RAPCEval ) 


Dennis  A.  Ludwig 
MSEE  Program 
UPDATE 


TECHNICAL  SECTION:  NEURAL  NETWORK  CONT. _ _ _ 

_ J  1.  GATHER  TEST*  TRAINING  DATA 

^AKl^i  2.  CHOOSE  INITIAL  WEIGHTS 


Designing  a 
Backpropagation 
Network 


ADJUST 

WEIGHTS 


TRAIN 

SYSTEM 

i 

_ 1 

N  Training  N 
N.  Completed 


ADJUST 

ARCHITECTURE 


TEST 

SYSTEM 


RF  RECEIVER  AND  PROCESSING  a.  Ludwig 

(al  LU  jp)  CONCEPTS  EVALUATION  PROGRAM  updat”" 

(RAPCEval)  ' - 


RF  RECEIVER  AND  PROCESSING 
CONCEPTS  EVALUATION  PROGRAM 
(RAPCEval) 


Dennis  A.  Ludwig 
MSEE  Program 
UPDATE 


TECHNICAL  SECTION:  NEURAL  NETWORK  CONT. 


Data  Collection  for  Training 

*  Set  model  in  pose  volume  and  measure  sampled 
points  from  focal  point. 

★  Re-orient  the  model  and  measure  sampled  points 
from  focal  point. 


★  Determine  the  new  pose  given  those  A-parameters 
(outputs)  and  the  resultant  minimum  distance 
vectors  (inputs) 


RF  RECEIVER  AND  PROCESSING 
(sf  rn  ]l)  CONCEPTS  EVALUATION  PROGRAM 

( RAPCEval ) 

TEC3iNICAL  SECTION:  NEURAL  NETWORK  CONT. 

Training  and  Testing  the  Network 


★  Data  is  divided  into  training  and  testing  sets 

★  Training  occurs  using  collected  data 

★  Trained  NN  is  then  evaluated 

★  Adjustments  made  if  results  are  not  acceptable 


Dennis  A.  Ludwig 
MSEE  Program 
UPDATE 


RF  RECEIVER  AND  PROCESSING 
CONCEPTS  EVALUATION  PROGRAM 
( RAPCEval ) 

PROPOSED  SCHEDULE 


Dennis  A.  Ludwig 
MSEE  Program 
UPDATE 


★  Quarter  1  r'An 

.  Develop  expertise  in  data  modeling  and  BRL-CAD 

software  (Didn't  use)  , 

.  Develop  expertise  in  neural  network  design 
package  (Didn't  use,  wrote  my  own) 

.  Develop  initial  solution  strategy  for  proposed 
problem 

★  Quarter  2  (Current  Semester) 

.  Develop  and  test  backpropagation  network 
.  Document  results 


I  could  get  the  model  to  move  to  the  image  by 
minimizing  an  angle  between  the  light  source 
ray  tracing  the  minimum  and  maximum 
points  of  the  two  images.  The  model  point  is 
well  defined,  but  the  Z  coordinate  will  be 
missing  from  the  target  image. 


Look  at  rotating  an  enclosed  area  to  get  Z 
information. 
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5.1.9  Presentation  by  Tracy  Tillman 

The  student  briefing  presented  by  Tracy  Tillman  at  this  meeting  is  reproduced  on  the  next  30 
pages. 
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Research  Topic: 

HARDWARE  IMPLEMENTATION  for  ADVANCED  PULSE  PROCESSING 


Quantify  Speed  Achieved  Through  Hardware  impien 

Establish  Hardware  Requirements 

Develop  Specification  for  liardware  Implementation 


..  Establish  Hardware  Requirements  .. 

..  Develop  Hardware  Specification . . 

Write  Project  Report . 

.  Submit  Project  Report 
Project  Complete  . . 
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SYMMETRICAL  ABOUT  1.5 

ODD  NUMBER  OF  ENTRIES  IN  BRANCH 
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Table  Entries 
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★  91  EQUIVALENT  FP  ADDITIONS 
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★  59  EQUIVALENT  FP  ADDITIONS 
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.  NUMBER  OF  OCCURRENCES  OF  EACH  TABLE  ENTRY 
OUTPUT:  VALUES  ABOVE  THRESHOLD 
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•  166.7  -  213.3 
.  126.7  -  246.7 
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•  REAL  TIME  EXPANSION  OF  TABLE 

•  LARGE  COUNTDOWN/SMALL  PRPs 

•  PULSE  TRAIN  ANALYSIS  WITH  LARGE  TABLES 


5.2  Steering  Committee  Meeting,  August  1998 
5.2.1  Meeting  Minutes 


EW  Receiver  and  Processing  Concepts  Evaluation  Program 
Program  Research  Standards  Committee  Meeting 
Minutes  August  12, 1998 


A  meeting  of  the  RAPCEval  Advisory  Committee  was  held  in  the  Mercer  Engineering  Research 
Center  Auditorium  at  1:00  p.m.,  August  12,  1998. 

The  meeting  consisted  primarily  of  an  invited  lecture  by  James  P.  Stephens,  Sr..  Dr.  Stevens 
works  in  the  RF  Sensor  Technology  Division,  Electronic  Warfare  Technology  Branch,  of  the  Air 
Force  Research  Laboratory  in  Dayton,  Ohio.  The  lecture  focused  on  the  topic  of  RF  sensor 
technologies  in  the  area  of  C3  Electronic  Combat. 

The  lecture  had  a  communication  countermeasure  emphasis,  focusing  on  advanced  digital  signal 
processing  for  EW.  It  began  with  an  organizational  background  and  objectives,  C2W  trends,  and 
vision  from  the  perspective  of  SNRW. 

A  brief  tutorial  was  provided  to  explain  the  motivation  of  the  use  of  time-frequency  distributions, 
wavelets,  and  cyclostationary  signal  processing.  Current  reseairch  activities  of  SNRW  in  all  areas 
related  to  C3  Electronic  Combat  were  introduced.  These  include  such  areas  as  featureless 
waveform  countermeasures,  time-frequency  analysis,  frequency-hopping  prediction,  GPS  ECM 
and  ECCM,  and  covert  waveform  detection  and  Jamming.  In-house  activities  include  time- 
difference  of  arrival  using  cyclostationarity,  frequency  hopping  synthesizer  development, 
spectral  correlation  analysis,  and  transform  domain  waveform  experimentation. 
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5.2.2  Attendance  Roster 


Bahl 


Bass 


Bass 


Black 


Bryant 


Campbell  Mark 


Erik  IMERC 


Charles  MERC 


Tom  MERC 


Joseph  K.  MERC 


Peter  MERC 


Carter 


DePuy 


Finnigan 


Guest 


Harris 


Hillhouse 


Jones 


McFarland 


Messer 


Miller 


Mitchell 


Moody 


Moody 


O'Dowd 


Oliver 


Pequignot 


Potts 


Glenn  IMERC  _ 


George  MERC  _ 


Skip  MERC _ 


Robert  MERC 


Rhondia  iMERC 


Ed 


Houston 


Joseph  F.  WR-ALC/LNIE 


Mark  MERC 


Chris  MERC 


Bob 


David 


Ray 


Douglas 


Douglas 


Bill 


R.P 


Nick 


Marvin 


IMERC 


WR-ALC/LNEV 


MERC _ 


SAIC 


SAIC _ 


Ball 


IWR-ALC/LNERT 


AFRL/SNRD 


AFRL/SNRW 


Smith 


Smith 


Thorstens  Gary 


Tillman  Tracy 


Wimpy  [Dale 


Jermaine  IMERC 


Jerome  WR-ALC/LNEV 


Ball 


MERC 


Phpne 

E-mail  V.,  -'  V 

912-953-6800 

ebahl@merc.mercer.edu 

912-953-6800 

cbass@merc.mercer.edu 

912-953-6800 

tbass@merc.mercer.edu 

912-953-6800 

jblack@merc.mercer.edu 

912-953-6800 

pbryant@merc.  mercer.edu 

912-926-5673 

mcamp- 

bel@avionics2.robins.af.mil 

912-953-6800 

gcarter@merc.mercer.edu 

912-953-6800 

gdepuy@merc.mercer.edu 

912-953-6800 

sfinnigan@merc.mercer.edu 

912-953-6800 

rguest@merc.mercer.edu 

912-953-6800 

rharris@merc.mercer.edu 

912-922-8106 

hillhouse@itc-2.com 

912-926-5458 

jhjones@hom.net 

912-926-4482 

kelley.joe@ec.robins.af.mil 

912-953-6800 

mmay@merc.mercer.edu 

912-953-6800 

cmcfarland@merc.mercer. 

edu 

912-953-6800 

bmesser@merc.mercer.edu 

912-926-2994 

miller@ec.robins.af.mil 

912-953-6800 

rmitchell@merc.mercer.edu 

912-918-2911 

moodyd@saic.com 

912-918-2911 

moddyd@saic.com 

912-922-4363 

bodowd@f1 5aisf.robins.af.mil 

912-926-2588 

oliver@ec.robins.af.m.il 

937-255-6127 

pequignot@sensors.af.mil 

937-255-5579 

pottsmn@sensors.wpafb.af. 

mil 

912-953-6800 

jsmith@merc.mercer.edu 

912-926-2994 

smithjm@ec.robins.af.mil 

912-922-4363 

gthorste@ba  11 .  com 

912-953-6800 

ttillman@merc.mercer.edu 

912-926-9529 

dwimpy@f1 5aisf.robins.af.mil 

5.2.3  Overview  of  the  Program  (Dr.  Bass) 

The  Overview  Briefing  of  the  RAPCEval  Program  as  presented  at  this  meeting  is  reproduced  on 
the  next  1 1  pages. 


149 


00 

o^ 

o^ 


^  ^  & 
O  w 

^  ^ 
u  E 


g|o 

>  w  w 

w  g  ^ 
u  H  Pm 

E  s  ^ 

^  o 
u 


o 

z 

I- 

UJ 

lU 

s 

"(S 

> 

LU 

O 

0. 

< 

0^ 


150 


August  12, 1998 
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committee  to  be  useful  to  me  Air  t  orce  ana  is 
certified  by  the  University  as  having  academic 
merit 
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available  from  MERC 
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of  Radar  Signals”,  available  from  MERC  on  need  to  know 

★  Kirk  Wright,  '‘Object  Oriented  Modeling  of  the  AN/ALQ- 
172”,  classified,  draft  in  edit  at  AFRL 
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C3  Electronics  Combat 


w 

a>  ^ 
O  C 

<u 


cd  o 

<u  ,o 

GO 

(D 

J-(  C/D 

o  2 

:S 

t4-l  E”* 

o  ^ 

4)  C/J 

•  ’-H 

E  ^ 

a  w 

o  ^ 

^  in 

cd  X2 


cd  2 

I  .£P 

c« 

,S  •- 

^  <=o 

QJ  • 

■C 

X)  ”0 

(D  <1^ 
^  O 

H  ^ 

d  -o 

O  as 

c 

•  *— <  ^ 

>  C) 

S.s 

GO 

bD  ;3 
O  o 

o  <S 

c  ^ 

H  CO 
O 

Cd 

H  ^ 
^  & 
CO  S 

a  ^ 

<D  CO 

GO  W 


(N  O 

c«  ^ 
(U  -o 

•  ^  > 

o  2 

•%  ^ 

^  CO 

O  •  '^ 


*0  , 
S  C  *=! 

Cu  •’— < 


•g  3 

1 

P  (i> 

bO  ^ 

X) 


oS 

c  ps^ 


®  B 

bO  W 

c  H 
A 

•g  = 

^  c/5 
'(— » 
CO  Cd 

H  I 

o 

U 


11 

C^-l  <u 

o  c 

^  bS 

c  8 

DQ  O 

bO  o 

^  C 
O  13 

C  P 

o  S 
w  o 
H  o 


os  CXJ 
N 

•p  O 
oS  <D 
bO  >  •. 

o  ^ 

■5  ^ 

•:r  00 

^  <D 
Ol  Oh 

.s  oT- 

bO  .£2  c 
(U  '►•' 

c 


i2  ^ 

p/  03 
<D  '7'  'O 

o  2 

tw  tw  -tb 

g  .3a  g 

.2  .ts  o 

3  :S 

-  2  2 
cb  JD 

.22  I 

^  2  T3 

>.  Cd  (U 

g  S 

p  2 

cr  c  -tt 

<U  D  c 

-V-t 

^  =  £ 
2  §  is 

.S  ts 

H  -2 

o  .  £ 


1:p  c 

W  .2 

U  CO 

c2  > 

bO  "2 

^  s 

*55 

s  4 

8  « 

2  p 


I 

(U 

£  '2 
.a  c 

cd 

-s 

CO 

2  ^ 

S3  U 

^  H 

p  w 

£  TS 


v-<  rj  — 

^  ^  8 

§  U  .2 

g  GO 
O  Pj  o 
O  cb 


w  O 
S£3  Jh 

a>  ^  C 


D  bb  o 

^  .£  u 

CW  o 
lU  CO  .P 

.£3  ^  C 
■^2  0 
G— I  P  S-i 

OJ-H  ^-‘ 

Cu  o 
C  h-h  a> 

.  2  g  W 

Cd  W)  rn 

o 

o 

£  I  ^ 

^  o  B 
■5  ■•5  i2 
c3  -g  sa 

cd  o  CO 
(U  O  as 


o  >— I 

Oh  bO 


^  Oh 
W  Oh 

cb  o 

00  h£h 


2P  •£ 

.s  £ 


C 

2 

(U  O 
Sh-2 

(U  o 
<1> 

-  h2 

CO  'O 

• 

CO  C 

£ 

g  ^ 

cb 


ITJ  S  ^ 

Qj  ;-i  .hh 

y  fP  > 

22  w  di 

1  ^ 

^  C  o 
bO  '§  '£ 

c  g 

'^1  ^ 

Oh'P  -p 

p  g  cb 

,0  Cd 
O  ^  n-t 

S  C  ^ 

2  £  'P 

^  v-i  si 

ST  <i> 

2  'p  ■£ 

c  .a 

^  to 

r^  •  f-H 

•C  'S  *30 

O  Cd  (D 
•  ^  cl  •  ^ 

tb  P  ^ 

■K  a  (u 

h2  .2  -g 

o  •+3  2 

.a 

o  o  a 

bO  £  • 

P!  9  ^ 

.-  o  OS 

p  13  .g 

>-|  .a 

^  O  H 
.a  (u 
H  P 

£  “  p 


.a  aJ 


H  P'  P 

p  “  .2 

o  P  .2 

^  p  c 

^  <-H  »-H 

o  a  p 
p  p,  P 
(U  O  -P 

Lh  Vh 

cu  CD  a> 

fe  >  Oh 

^  B  X 

’Z3  'O  <o 


159 


current  activities  of  SNRW  and  should  be  of  interest  to  technical  managers, 
researchers,  faculty,  and  engineering  students. 
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in  a  three-state  area  with  prior  assignments  in  Detroit  and  Norfolk.  Since  Jan  1982,  Mr.  Stephens  has  been 
employed  part-time  as  an  instructor  in  the  Electrical  Engineering  Technology  curriculum,  Associate  Degree 
program,  at  Miami  University  where  he  teaches  a  variety  of  courses  in  electronic  and  microprocessor 
technology.  He  has  authored  or  co-authored  more  than  15  publications. 


5.2.4  Presentation  by  Jim  Stephens,  AFRLVSNRW,  Dayton 


A  guest  lecture  presented  by  James  P.  Stephens,  Sr.  at  this  meeting  was  described  by  a  synopsis 
in  the  overview  of  the  previous  section.  Copies  of  this  presentation  are  available  on  request  from 
the  RAPCEval  Program  Manager  at  Mercer  Engineering  Research  Center. 
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5.3  Steering  Committee,  November  1998 
5.3.1  Meeting  Minutes 


EW  Receiver  and  Processing  Concepts  Evaluation  Program 
Program  Research  Standards  Committee  Meeting 

Minutes  9  November,  1998 _ 


The  RAPCEval  Steering  Committee  held  a  meeting  at  MERC  on  9  November,  1998  at  1 :30 
p.m.  Dr.  David  Barwick,  MERC  Executive  Director,  kicked  off  the  meeting  with  a  welcome. 

Dr.  Tom  Bass,  MERC  RAPCEval  Program  Director,  then  presented  a  “State  of 
RAPCEval”  slide  overview  presentation. 

Dr.  Bass  then  introduced  the  three  students  who  were  to  give  presentations  on  their  research 
topics.  Those  students  and  their  topics  were: 

♦  Mr  Houston  Jones  -  “Evaluation  of  Reed-Solomon  Codes  for  CDMA  Systems" 

♦  Mr  Mark  Napier  -  “Civil  IFF  Reed-Solomon  Code  Application” 

♦  Mr.  Randy  Ford  —  “Passive  Location  via  Evolutionary  Genetic  Algorithms” 

The  next  agenda  item  was  a  period  of  discussion  and  suggestions  about  the  projects.  Ches 
Rehberg  (during  discussions  and  via  a  later  e-mail)  pointed  out  the  need  for  Air  Force  personnel 
adept  in  the  topics  being  presented  to  attend  the  meetings.  For  example,  two  of  the  topics  were 
on  communications  (robustivity  via  Reed-Solomon  code  implementations).  Ches  suggested  such 
participants  for  future  meetings. 

Phil  Oliver  stated  that  data  links  are  a  big  topic  in  the  F-15  commumty.  He  suggested  that 
these  communication  robustivity  research  topics  may  have  an  application  in  that  area. 

Several  attendees  remarked  that  the  student  presentations  appeared  higher  in  quality  on  this 
occasion.  Some  questions  were  raised  about  the  applicability  of  the  communication  topics  to  the 
original  RAPCEval  charge. 

The  meeting  was  then  adjourned. 
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5.3.2  Agenda 


Meeting  Agenda,  continued,,. 


Schedule 

Greetings 

Dr  David  Barwick 

1:30  PM  ^1:35  PM 

Meeting  Overview 

Dr  Tom  Bass 

1:35  PM -1:45  PM 

Student  Proposals 

;Sfodents: 

Evaluation  of  Reed- 
Solomon  Codes  for 
CDMA  Systems 

Mr  Houston  Jones 

1:45  PM -2:15:  PM 

Civil  IFF  Reed-  1 

Solomon  Code  Ap-  : 
plication 

Mr  Mark  Napier 

2:15  PM -2:45  PM 

Passive  Location  via 
Evolutipiiary  Ge¬ 
netic  Algorithms 

Mr.  Randy  Ford 

2:45  PM -3:15:  PM 

Discussions  &  New 
Business 

Dr.  Tom  Bass 

3:15  PM -3:45  PM 

A(^oum 

4:15  PM 
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5.3.3  Attendance  Roster 


The  attendees  at  this  meeting  are  listed  here; 


Name 

Organization 

Phone 

Email  Address 

I 

Tom  Bass 

MERC 

912-953-6800 

bass_wt@mercer.edu 

H 

Dave  Barwick 

MERC 

912-953-6800 

barwick_dt@mercer.edu 

■ 

Ches  Rehberg 

WR-ALC 

/LNEX 

912-936-4525 

rehberg@ec.robins.af.mil 

Kl 

Charles  Bass 

MERC 

912-953-6800 

bass_cd@mercer.edu 

■ 

RP  (Phil) 

Oliver 

WR-ALC 

/LNERT 

912-926-2588 

oliver@ec.robins.af.mil 

■ 

Steve  Strawn 

WR-ALC 

/LNERR 

912-926-6435 

strawn@ec.robins.af.mil 

■ 

Tony  White 

AFRL 

/SWRP 

937-255-6127, 
Ext.  4236 

whiteaw@sensors.wpafb.afmil 

8 

Nicholas 

Pequignot 

AFRL 

/SNRP 

937-255-6127, 
Ext.  4235 

pequigna@sensors.wpafb.af.mil 

9 

Mark  Napier 

Scientific  Atlanta 

770-903-6980 

mark.napier@subasic.sciatl.com 

■ 

Randy  Ford 

WR-ALC 

/LYSBC 

912-926-0423 

fordra@avionics.robins.af.mil 

■ 

Houston  Jones 

WR-ALC 

/QLYM 

912-926-5458 

jhjones@hom.net 

i 

Behnam 

Kamali 

Mercer 

912-752-2415 

kamali_b@mercer.edu 

13 

Paul  McNeil 

Mercer 

912-752-2185 

macneil_pe@mercer.edu 

14 

Joseph  K. 

Black 

MERC 

912-953-6800 

jblack@merc.mercer.edu 

15 

Ray  Mitchell 

MERC 

912-953-6800 

rmitchell@merc.mercer.edu 

16 

Skip  Finnigan 

MERC 

912-953-6800 

sfinnigan@merc.mercer.edu 
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5.3.4  Overview  of  the  Program  (Dr.  Bass) 

The  Overview  Briefing  of  the  RAPCEval  Program  as  presented  at  this  meeting  is  reproduced  on 
the  next  7  pages. 
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to  be  useful  to  the  Air  Force,  as  well  as  having  merit 
for  graduate  research 
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*  MERCER  UNIVERSITY:  Dr.  Aaron  Collins,  Dr. 
Benham  Kamali,  Dr.  Paul  MacNeil 
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Solomon  Codes”  to  be  included  in  forthcoming  unclassified 
report 

★  Mark  Campbell,  “Auto-Regressive  Spectral  Analysis  - 
EW  Applications”,  unclassified  report  available  from  MERC 
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5.3.5  Presentation  By  Mark  Napier 

The  student  briefing  presented  by  Mark  Napier  at  this  meeting  is  reproduced  on  the  next  1 5 
pages. 
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APPLICATION  OF  REED-SOLOMON  ENCODING  TO  IMPROVE  PROPOSED 
COLLISION  AVOIDANCE  SYSTEM  BASED  ON  CIVILIAN  ATCRBS 


The  new  message  is  112  bits  long,  40  of  which  have  been  reserved  for  FEC 
coding.  The  proposed  FEC  scheme  is  a  5  bit  t=4  Reed-Solomon  code, 
RS(31,23).  The  proposed  work  is  to  analyze  this  scheme  to  determine  the 
reliability  improvements  to  be  realized  from  the  RS  code.  A  decoder  will  be 
designed,  implemented  and  tested  in  Verilog. 


•  Similar  encoding  could  be  used  in  Military  IFF 

.  Military  aircraft  could  participate  in  the  civilian 
system  without  revealing  aircraft  ID 
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Design  the  Decoder  in  Verilog 
Verify  with  Verilog  Testbench 
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★  Model  the  System  with  RS  FEC  coding 

•  Model  without  Erasure  Information 

•  Model  With  Erasure  Information 
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Develop  and  test  decoder  for  RS(31,23)  code. 
Document  results  and  submit  project  report. 
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Position  Modulation)  scheme  with  a  1  Mbit/s  rate.  A  “1”  is  defined  to  be  a  0.5us  burst 
followed  by  0.5us  of  off  time.  A  “0”  is  defined  to  be  0.5us  of  off  time  followed  by  a  0,5us 
burst.  The  message  is  proceeded  by  a  8us  sync  pulse.  Either  56  (single  length)  or  1 12  (double 
length)  bits  of  data  follow. 
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Mode  C  Packet  -  Identical  to  Mode  A  packet.  Altitude  encoded  on  10  bits  of  the  digit  values. 
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parity  bit  for  message. 


Schedule  graduation  (campus) 
Project  committee  defense 
Submit  project  report  to  university 


5.3.6  Presentation  by  Houston  Jones 

The  student  briefing  presented  by  Houston  Jones  at  this  meeting  is  reproduced  on  the  next  1 5 
pages. 
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SOLOMON  CODES  IN  WIRELESS  CODE  DIVISION  MULTIPLE  ACCESS 
SYSTEMS 
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Digital  cellular  communications  are  important  and  valuable  to 
military  as  well  as  civilian  users.  Increases  the  knowledge  base 
in  digital  mobile  communications 
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communications  systems 

Determine  the  feasibility  of  using  Reed-Solomon 
(RS)  codes  to  reduce  system  complexity  in 
CDMA  mobile  communications  systems 
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error  rate 

Optimistic:  Improved  performance  and  reduced 
complexity 
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Performance  --  RS  <  CC+BI; 

RS  =  CC  +  BI  or  RS  >  CC  +  BI 
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Study  to  determine  the  feasibility  of  using  RS 
code  families  for  single  burst  error  correction  in 
Rayleigh  fading  channel 
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Simulate  CDMA  system  with  convolutiona 
encoder  and  block  interleaver 
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Error  correction  capability 
Frame  error  rate;  bit  error  rate 
Delay 

System  complexity 
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In  1993  a  cellular  system  based  on  code  division  multiple  access 
(CDMA)  was  developed.  This  system  uses  direct  sequence  spread 
spectrum  techniques  requiring  less  power  and  allowing  for  more 
•■•userS"-than-the-other"designS'. . 
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Forward  CDMA  channel  modulation  process. 


5.3.7  Presentation  by  Randy  Ford 


The  student  briefing  presented  by  Randy  Ford  at  this  meeting  is  reproduced  on  the  next  47  pages. 
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esearcn  lopic: 

COMPARISION  OF  DIFFERENTIAL  EVOL  UTION  TO  THE  SIMPLEX  METHOD  IN 
OPTIMIZA  TION  DURING  PASSIVE  EMITTER  LOCATION. 
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Differential  Evolution,  developed  by  Kenneth  Price  and  Rainer  Storn,  is  a 
minimization  technique  for  real-valued  problems  which  is  based  on  genetic  algorithms. 
This  project  will  compare  the  performance  of  this  method,  with  the  performance  of  the 
Nelder-Mead  simplex  method.  Both  methods  will  be  implemented  in  C++  and  called 
from  an  existing  MATLAB  implementation. 
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Test  both  implementations 

•  Identical  Platforms  -  Pentium  PC 

•  Determine  best  method  based  on  speed  of  detection. 
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x,y  plane. 
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★  Passive  Emitter  Location  discussion  is  based  on  An 
Efficient  Method  of  Passive  Emitter  Location”  by  Klaus 
Becker  in  the  October  1992  IEEE  Transactions  on 
Aerospace  and  Electronic  Systems. 
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Reflection 

and 

Expansion 


I 

I 

I 

I 

i 

til 

i 

I 

I 

I 

I 

I 


I 

o 

I 


Et^ 

S 

§ 

Co 

§ 

I 

I 


o 


H 

< 

H 

Iz; 

w 

C/) 

S 

flH 

H-l 

< 

o 


u 

w 

H 


o 

a. 

jo 

a 


c/3 

C 

o 

»  F-^ 

CO 

n 

<D 

6 


c3 


u 

DD 

CD 

•  fH 

W 

X! 

1 

Jh 

s 

•  f-H 

o 

CO 

rS 

03 

CIh 

o 

D 

C3 

43 

CO 

31 

D 

o 

O 

CCJ 

■ 

H 

• 

218 


Low  All  Dimensions 
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Differential  Evolution  is  an  Evolutionary  Strategy 


a  00 

OTs 

o  OX) 

^  2  ^ 

-a 

c«  w  S 

(1>  C/!)  W 

I  W  o 


Xfl 

H 

pH 

R2 


5  psj  « 

^  ^ 

S  u 

^  O  Ph 

y  p  ^ 

Iz;  nJ 

O  ^ 

A  > 

H  [id 


1  ^f" 


i  ys 

I# 


1^ 


^2 


^  s 

gg 


3l 

&  K 

fcl 


x-<  ^ 

So 

S 


s  - 

Oh  <D 

•  •  bO 

c/)  ^03 

(D  .S  r^ 

•’-'  tii  Cu 

t>o  s 

g  C/1  .2 

^  ^  s 

H  c^  ^ 

3  >  O 

§  ^  'S 

O  ^  -f-H 

8  S 

—  c  a 

o  w  ^ 


(D  ^ 

g  T3 

2^  c^o 

'S 

GO  pH 

S  > 

*,j2  ^ 

o  ^ 


^  -2 

S 

’-p 

x/1  c/3 


C  ^ 

C  c3 

(D  (D 

4-»  P 

<D  c/3 

P  P 

^  P 

"  o 


P  (U 

c/3  C/3 

;d  P 


222 


This  makes  incremental  search  easier. 


in 


Ph 

h-i 


o 


o 

> 

PQ 

•  ^pH 

'S 

(D 

cS 


c/j 

O 


<D 

> 

I 

<D 

Vh 


a 

o 

►  F-^ 

'S 

N 

*  F^ 

►  f— H 

4-> 

►  t-H 

a 


c/3 

a> 

Id 

> 

;-i 

<D 

4-* 

<D 


a 

o 

c/3 

cd 

(D 


o 

cd 

(D 

4h 

O 

'tNk 

l>3 

(D 

^  . 

■cd  ^ 


2 

> 


Cd 

g 


T3  ^ 

S 

cd 

o  C 

O  •  F— 4 
'■P  CO 

cd  QJ 

3  3 

Ph  ^ 
O  > 

D 

»— H  C/3 
Cd  CD 

:s  ■s 

•S  j> 

P  <D  O 


cd 

c3 

P. 

<u 

P 

o 

CO 


CD 


2  S 

S  ^ 
^  S 

O  g 


223 


Each  member  of  the  current  population  takes  its  turn  as  the 
target  vector  for  mutation,  crossover,  and  selection 
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Mutation: 

■  Determine  the  difference  between  two  randomly  chosen 
population  members,  x^  -  x^ 

>  Multiply  by  the  scaling  factor  F.  F*(  x^  -  x^,) 
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parameter. 


i-i 

O  <50  ^ 

^  B  ^ 
^  -g 

CO  W  S 
(U  C/D  ^ 

I  W  o 
^  22  ^ 


Tfi 

'H 

8^ 

w 

^  ^  5 

<;  oj  > 

^  w 

z  u 

■<  O  Pm 

Z  1-5 

o 

Oi  >; 

H  (i5 


as 

^  Ci 

l§ 

ii 

cq 

tS 

|g 

^5 

as 

fcR 

5 

Ci  N 

®  c 
So 
aa  ^ 


5^  • 

c3  ^ 

Qm  O 

B  ^ 

o  o 

s 

a>  ^2 

^  ?-H 

S  ^ 

S  rO 

fe  B 

^  § 
^  S 

^  B 

03  S 
.  ^ 

I  § 


'i  S 

2  ^ 
,13 

<D 

^  S 

^  P 

HH 


;-H 

<U 

i 


226 


:rom  x^  . 

rhe  last  parameter  always  comes  from  x^,’  in  order  to  make  sure 
hat  the  trial  vector  differs  from  the  target. 
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Trials  continue  for  a  pre-set  number  of  iterations  or  until  the 
cost  of  the  best  vector  is  less  than  or  equal  to  a  pre-set  value. 
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Next  Generation 
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Calls  external  C++  dlls  called  MEX  files  for 
minimization 

Two  versions  -  nmloc.m  and  deloc.m 
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Version  deloc.m  calls  C++  file  dec95.dll  to  minimize 
using  DE. 

Both  read  in  data  from  same  external  file  data.m 


scaling  factor,  crossover  constant,  and  number  of 
iterations  from  the  user. 
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data.m  floatf.exe  nmc275.dll 
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distances, 
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Searches  the  error  space  by  calling  de_minimizer. 
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Use  the  same  data  while  adjusting  the  parameters  of  the 
DE  simulator  to  try  to  improve  its  performance. 
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and  more  accurate  in  only  5  of  the  12  test  sets. 
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Repeating  tests  with  the  same  input  data  results  in  the 
same  AED. 
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Further  study  of  parameter  combinations  may  make  DE 
even  more  accurate. 

Further  study  of  alternate  crossover  techniques  may 
make  DE  faster  and  more  accurate 
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5.4  Thesis  Report  by  Ron  Brinkley 

A  version  with  standard  report  style  has  been  generated  from  the  Master’s  thesis  of  Ron 
Brinkley.  This  is  reproduced  here  on  the  next  83  pages. 
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1  Summary 

Communication  and  storage  systems,  whether  they  are  magnetic  or  optical  recording  de¬ 
vices,  satellite  or  cellular  mobile  communication,  or  even  military  communication  sys¬ 
tems  in  hostile  jamming  environments,  are  faced  with  a  need  to  correct  burst  errors.  It  is 
the  intent  of  this  research  work  to  explore  the  efficiency  of  powerful  burst  error- 
correcting  Reed-Solomon  (RS)  codes.  First,  an  explanation  as  to  the  cause  of  burst  er¬ 
rors,  and  a  definition  of  burst  errors  is  presented.  This  explanation  is  followed  by  an  in¬ 
troduction  to  the  finite  field  mathematics  necessary  in  understanding  RS  codes.  An  in¬ 
troduction  to  RS  codes  and  their  ability  to  combat  burst  errors  is  given.  Subsequently  an 
overview  of  the  most  widely  used  interleaving  techniques  is  presented  to  complete  the 
preliminary  portion  of  this  thesis. 

The  guaranteed  biurst  error-correcting  capability  of  RS  codes  is  first  explored  for  sin¬ 
gle  and  multiple  bmrst  errors.  This  burst  error-correcting  capability  is  then  used  to  define 
the  burst  error-correcting  efficiency  of  selected  RS  codes.  The  idea  of  guaranteed  burst 
error  correction  and  burst  error  correction  efficiency  is  then  expanded  to  include  block 
and  convolutional  interleaving  techniques.  Many  useful  tables  and  graphs  of  guaranteed 
burst  error-correcting  capabilities  and  efficiencies  have  been  generated  for  the  code  de¬ 
signer. 
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2  Introduction 


Why  the  concern  for  burst  error  correction?  Other  than  deep  space  communications, 
an  overwhelming  majority  of  practical  communication  systems  carmot  be  modeled  as 
simple  Additive  White  Gaussian  Noise  (AWGN)  channels.  In  recording  systems, 
whether  they  are  magnetic  or  optical,  defects  in  the  recording  material  and  recording 
mechanisms  cause  systematic  errors.  Wireless  communication  systems  with  fading  ef¬ 
fects  are  also  faced  with  the  occurrence  of  non-random  errors.  Jamming  environments, 
especially  intentional  jamming  as  in  military  applications,  are  confronted  with  the  effects 
of  non-random  errors.  Communication  channels  that  have  non-random  characteristics  are 
considered  to  contain  memory  and  tend  to  produce  what  is  known  as  burst  errors.  In  fact, 
most  operational  communication  system  applications  exhibit  a  combination  of  random 
and  non-random  characteristics.  These  channels  are  termed  diffuse  or  compound  chan¬ 
nels  [Swe]  and  cannot  be  modeled  as  AWGN  channels.  Therefore,  in  answering  the 
question  why  burst  error  correction,  most  real  world  communication  systems  need  codes 
that  can  combat  the  common  occurrence  of  burst  errors. 
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2.1  Burst  Errors 


The  purpose  of  this  section  is  to  provide  the  reader  with  a  better  understanding  of  the 
errors  associated  with  digital  communication  systems.  First  a  simple  example  and  defi¬ 
nition  of  errors  that  occur  in  binary  codewords  is  given.  A  discussion  of  the  two  different 
types  of  possible  errors  follows.  Then  the  definition  of  errors  is  expanded  to  include  the 
word  oriented  RS  codes.  The  section  is  concluded  with  brief  examples  of  the  classic  er¬ 
ror  causing  bursty  channel. 

A  received  encoded  block  of  data  can  be  interpreted  as  the  linear  combination  of  a 
valid  codeword  and  an  error  pattern.  This  relationship  is  shown  as  R  =  C  +  E  where  R  is 
equal  to  the  received  block  of  data,  C  is  a  valid  codeword  and  E  is  an  error  pattern  and  the 
summation  implies  binary  vector  addition.  In  the  binary  case  a  ‘  1  ’  in  the  error  pattern 
represents  the  occurrence  of  an  error  in  that  position.  As  an  example,  consider  a  commu¬ 
nication  system  where  the  sixteen  bit  block  of  received  data  is  R  =  0010001 1 10010000. 
Then  suppose  the  original  transmitted  data  was  the  valid  codeword  C  = 
0001010010001000.  Therefore,  because  of  errors  in  the  channel,  an  error  pattern  E  = 
001101 110001 1000  is  produced.  It  can  be  observed  that  the  modulo  two  addition  of  R 
and  E  will  correct  the  errors  in  R  and  return  the  intended  block  of  data  C. 

Errors  occur  in  communication  systems  in  two  different  forms,  random  and  burst  er¬ 
rors.  An  example  of  a  random  error  is  given  in  Figure  2.0  where  n  bits  of  a  received  data 
stream  are  given  and  x  represents  a  bit  error.  There  are  five  random  errors  shown  in  Fig¬ 
ure  2.0  and  the  error  pattern  would  be  E  =  0001010100011000  with  weight  e  equal  to 
five.  Note  that  for  no  known  deterministic  reason,  two  of  the  random  errors  occurred  ad- 
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jacent  to  each  other.  Additive  white  Gaussian  noise  channels  are  random  error  channels. 
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Figure  2.0  Random  Errors 

Burst  errors  are  conventionally  defined  in  terms  of  their  length,  and  can  be  viewed  as  a 
cluster  of  individual  errors.  The  length  of  a  burst  error  is  equal  to  the  number  of  bits  that 
span  the  error  cluster.  Therefore,  a  burst  error  of  length  b  is  defined  as  the  number  of  bits 
in  the  error  cluster  starting  from  the  first  bit  in  error  (first  ‘  1’  in  the  error  pattern)  and 
ending  with  the  last  bit  in  error  (last  ‘  1’  in  the  error  pattern)  [Lin].  For  example  the  burst 
error  shown  in  Figure  2.1  is  of  length  b  =  1 1  bits  and  has  an  error  pattern  of  E  = 

001101 110001 1000  with  weight  e  equal  to  seven. 
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Figure  2.1  Conventional  Burst  Error 

More  than  one  burst  error  may  be  imposed  onto  a  codeword.  When  this  happens,  the 
codeword  is  said  to  contain  multiple  burst  errors.  For  example  in  Figure  2.1  above,  the 
codeword  may  be  considered  to  contain  three  burst  errors  of  which  two  are  of  weight  two 
and  one  burst  error  of  weight  three.  In  keeping  with  the  notation  used  in  [Ben],  a  code¬ 
word  may  contain  s  burst  errors  (s  >  1)  of  length  b  bits  (b  >  1)  and  each  burst  of  weight  e 
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(e>l). 


To  resolve  any  ambiguities  about  defining  burst  errors,  the  idea  of  a  guard  space  is 
now  explained.  Just  as  a  charmel  with  memory  produces  error  clusters,  the  channel  also 
produces  spaces  where  no  errors  occur.  These  clusters  of  error  free  bits  are  indicated  as 
sequences  of  all  zeros  in  the  error  pattern  E  and  are  called  guard  spaces.  Therefore,  a 
burst  error  b  is  defined  in  relation  to  the  surrounding  guard  space  g  [Gall].  For  example, 
in  Figure  2.2  the  knowledge  of  the  guard  spaces  dictates  the  definition  of  the  possible 
burst  errors.  The  bits  of  information  indicated  by  x’s  in  Figure  2.2  can  represent  burst 
errors  of  length  less  than  or  equal  to  b. 
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Figure  2.2  Burst  Errors  with  Guard  Spaces 

The  definition  of  burst  errors  will  now  be  expanded  for  the  case  of  RS  codes  that  are 
examined  in  this  report.  Reed-Solomon  codes  are  word-oriented  codes,  in  which  every 
codeword  is  composed  of  m  bit  symbols.  We  assume  that  the  random  error-correcting 
capability  of  a  RS  code  is  denoted  by  t.  If  any  one  of  the  m  bits  of  a  given  symbol  is  in 
error,  then  the  entire  symbol  is  considered  to  be  erroneous.  Therefore,  it  is  customary  to 
measure  random  errors  in  a  RS  code  according  to  the  number  of  symbols  in  error  and  not 
the  number  of  bits.  Figure  2.3,  a  representation  of  a  n  symbol  codeword,  indicates  two  (t 
=  2)  random  errors. 
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n  symbols 


There  are  several  systems  of  interest  that  operate  over  the  classic  bursty  channel  which 
use  burst  error-correcting  RS  codes.  Burst  errors  are  found  in  both  magnetic  and  optical 
recording  systems.  In  digital  mobile  communication  systems,  such  as  cellular  and  satel¬ 
lite  communication  systems,  burst  errors  are  a  major  concern.  And  in  military  communi¬ 
cations,  intentional  jamming  may  introduce  burst  errors  into  the  channel.  Burst  errors  can 
have  serious  adverse  effects  on  signal  reliability,  and  hence  the  need  for  burst  error-cor¬ 
recting  codes. 

With  increases  in  recording  speeds  and  data  storage  densities  for  magnetic  recording 
systems,  burst  errors  have  become  common  place.  Dust  and  defects  in  the  recording  ma¬ 
terial  and  the  read/write  heads  of  a  magnetic  disk  recording  system  are  usually  to  blame 
for  the  occurrence  of  bxarst  errors.  These  defects  have  been  known  to  produce  burst  errors 
of  hundreds  of  bits  in  length  [Lin,  chap  16].  One  example  of  the  use  of  RS  codes  in  mag¬ 
netic  recording  systems  is  the  IBM  3370  disk  storage  system  [Meej. 

Optical  recording  systems  have  used  RS  codes  extensively  for  the  correction  of  burst 
errors.  One  digital  recording  system  that  owes  its  existence  to  RS  codes  is  the  compact 
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disk  (CD)  recording  system.  With  the  reduction  in  size  of  the  recording  media  (a  com¬ 
pact  disk  has  an  outer  diameter  of  120mm  [Car]),  the  density  of  the  data  has  become 
enormous.  The  track  width  of  a  CD  is  roughly  0.6  pm  and  one  bit  of  audio  data  covers 
1 .0  pm^  of  disk  space  [Car].  In  comparison  to  the  LP  record  track,  the  CD  has  a  track 
ratio  of  approximately  60:1  [Bae].  That  is,  sixty  tracks  can  be  placed  on  a  CD  in  the 
same  amount  of  space  once  used  by  one  LP  track.  With  these  data  densities,  it  is  clear 
that  the  slightest  defect  or  scratch  may  cause  extremely  long  burst.  If  it  were  not  for  the 
error  control  scheme  known  as  Cross  Interleaved  Reed-Solomon  Code  (CIRC),  the  Com¬ 
pact  Disk  recording  system  could  not  deliver  nearly  the  quality  of  music  it  now  does. 

The  CIRC  error  correction  scheme  has  the  capability  to  correct  burst  errors  of  approxi¬ 
mately  4000  bits  in  length,  this  corresponds  to  roughly  2.5  mm  of  track  length  [Inun2]. 
For  a  complete  description  of  the  CIRC  error  control  method  as  applied  to  the  compact 
disk  refer  to  [Vri]  and  [Imml]. 

Since  the  onset  of  the  digital  communication  era,  mobile  communication  has  become  a 
fast  developing  area.  Popular  topics  in  the  mobile  cormmmication  arena  include  cellular 
telephone  systems,  personal  communication  systems  (PCS),  and  satellite  communication 
systems.  All  of  these  communication  systems  are  faced  with  channel  problems  such  as 
multipath  effects,  channel  fading  and  Doppler  shifts.  In  light  of  such  adverse  channel 
impairments,  burst  errors  are  very  common.  One  example  of  the  use  of  RS  codes  in  mo¬ 
bile  communications  is  in  the  cellular  digital  packet  data  (CDPD)  system  [Kami] 

[Kam4].  For  a  more  detailed  review  of  the  signal  characteristics  of  a  mobile  communi¬ 
cations  channel,  see  [Vau]. 
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Fading  in  a  mobile  communication  channel  can  be  statistically  modeled  as  a  Rayleigh 
distributed  process  [Wic3].  Fades  are  categorized  into  what  are  known  as  slow  and  fast 
fades.  The  size  of  the  burst  error  produced  is  directly  related  to  the  length  of  the  fade  im¬ 
posed  on  the  system.  Slow  fades  are  longer  in  duration  relative  to  the  symbol  rate  of  the 
communication  signal  being  affected.  Therefore,  a  slow  fade  will  produce  multiple  sym¬ 
bol  errors.  It  has  been  shown  that  mobile  communication  systems  v^dth  the  aid  of  error 
control  coding  and  interleaving  techniques  have  improved  performance  [Yue]. 

In  today’s  military  combat  environment,  the  ability  to  communicate  determines  win¬ 
ners  or  losers.  For  this  reason  electronic  warfare  (EW)  and  information  warfare  (IW)  are 
targets  of  many  discussions  and  much  research.  Of  particular  interest  is  the  survival  of  a 
communication  system  in  the  presence  of  intentional  jamming.  There  are  many  different 
types  of  jamming  techniques,  however  pulsed  jamming  has  shown  to  be  the  most  optimal 
for  creating  interference  [Tor].  Pulsed  jamming  concentrates  more  of  the  total  signal 
power  into  each  pulse  unlike  a  continuous  wave  jammer  whose  power  is  more  spread  out. 
If  the  pulse  width  of  the  jamming  signal  is  of  the  same  order  as  the  symbol  transmission 
rate,  then  approximately  one  symbol  will  be  affected  per  jamming  pulse.  If  the  pulse 
width  is  greater  than  the  symbol  transmission  rate,  then  multiple  symbols  can  be  affected 
per  jamming  pulse  causing  burst  errors.  As  illustrated  in  Figure  2.6  below,  jamming  sig¬ 
nal  J,(t)  has  pulse  width  approximately  equal  to  the  symbol  rate  of  the  communication 
signal  S(t).  Jamming  signal  J2(t)  has  a  pulse  width  approximately  equal  to  three  times 
the  symbol  rate  of  S(t)  thus  causing  burst  errors. 
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Figure  2.6  Pulsed  Jamming 

Techniques  such  as  spread  spectrum  communication  have  been  employed  with  good 
success  in  defending  against  intentional  jamming.  However,  it  can  be  shown  that  using 
spread  spectrum  with  the  addition  of  error  control  coding  and  interleaving  techniques,  the 
affect  of  pulsed  noise  jeimming  can  almost  be  eliminated  [Zie,  chap  1 1&12]  [Pur].  One 
example  of  a  military  communication  system  that  uses  RS  codes  is  the  Joint  Tactical  In¬ 
formation  Distribution  System  (JTIDS)  being  supported  by  all  branches  of  the  military 
and  NATO  [Zie,  chap  13]. 

With  the  ever  increasing  complexity  of  communication  systems,  the  powerful  burst 
error-correcting  RS  codes  will  no  doubt  see  more  applications  in  digital  recording  sys¬ 
tems,  mobile  communication  systems,  and  in  communication  environments  with  interfer¬ 
ence. 
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3  Finite  Field  Mathematics 


To  understand  Reed-Solomon  codes  an  introduction  to  Galois  (pronounced  ‘gal-wa’) 
field  algebra  is  presented.  To  learn  and  understand  Galois  fields  it  is  important  to  have  a 
foundation  in  the  laws  of  abstract  algebra  that  govern  Galois  field  theory.  However,  it  is 
by  no  means  the  intention  of  this  section  to  present  an  in-depth  coverage  of  abstract  alge¬ 
bra.  The  purpose  of  this  section  is  to  present  an  overview  of  basic  ideas  and  laws  needed 
in  understanding  Galois  field  algebra  as  applicable  to  RS  coding.  The  reader  is  encour¬ 
aged  to  review  any  of  the  following  references  for  a  more  detailed  analysis  of  abstract 
algebra  and  Galois  field  theory  [Cla],  [Lin],  [McEl],  [Petl],  and  [Wicl]. 

3.1  Groups 

A  group  (G)  is  a  set  in  which  a  binary  operation  ( •  )  has  been  defined.  Binary  in  the 
sense  that  the  operation  is  performed  on  two  elements.  The  resultant  of  the  binary  opera¬ 
tion  is  also  another  element  in  the  set.  This  unique  property  is  called  closure.  Three 
other  properties  are  required  in  order  for  a  set  to  form  a  group.  The  law  of  associativity 
must  be  upheld.  There  must  exist  an  identity  element  in  the  group.  Every  element  in  the 
group  must  have  a  unique  inverse.  Optionally,  if  a  group  satisfies  as  a  fifth  property  the 
law  of  commutativity,  the  group  is  said  to  be  commutative  and  called  an  Abelian  group. 
In  summary,  the  follow  properties  apply: 

a)  Closure:  V  a,b  e  G;  a  •  b  =  c;  where  c  €  G 

b)  Associativity:  (a  •  b)  •  c  =  a  *  (b  •  c) 

c)  Identity:  V  a  s  G  there  exist  e  €  G  such  that  a  *  e  =  e  •  a  =  a 
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d)  Inverse:  V  a  e  G  there  exist  a'*  e  G  such  that  a  •  a''  —  a’'  •  a  -  e 

e)  Commutativity:  V  a,b  €  G;  a  •  b  =  b  •  a  (optional) 

3.2  Rings 

A  ring  ( R )  is  a  set  of  elements  with  two  binary  operations  ( +  and  • )  defined.  There 
are  three  properties  that  must  hold  for  a  set  to  form  a  ring.  The  ring  must  form  an  Abe¬ 
lian  group  xmder  (  + )  with  the  additive  identity  element  indicated  by  (  0 ).  The  law  of 
associativitity  must  be  upheld  for  the  binary  operation  ( •  ).  Also,  the  binary  operation  of 
( • )  must  be  distributive  over  the  binary  operation  of  ( + ).  Additionally,  if  the  operation 
(  •  )  is  commutative,  then  the  ring  is  called  a  commutative  ring.  If  the  operation  ( • )  has 
an  identity  element  labeled  as  ( 1  ),  then  the  ring  is  said  to  be  a  ring  with  identity.  If  the 
ring  meets  both  of  the  additional  properties,  then  it  is  called  a  commutative  ring  with 
identity.  The  ring  is  an  important  building  block  for  error  control  coding.  It  is  the  ring 
that  allows  for  the  development  of  the  polynomial.  Summarizing,  the  following  proper¬ 
ties  apply: 

a)  Commutative  group  under  addition  (  + )  with  identity  element  (  0 ) 

b)  Associativity  imder  multiplication  ('):  V  a,b,c  e  R;  (a  •  b)  *  c  =  a  •  (b  •  c) 

c)  Multiplication  ( • )  distributes  over  addition  ( + ):  V  a,b,c  e  R; 
a  •  (b  +  c)  =  (a  •  b)  +  (a  •  c) 

d)  Commutativity  under  multiplication  (•)-'^^bsR;  a-b  =  b-  a 

e)  Multiplication  ( • )  has  an  identity  element  labeled  as  ( 1  ) 

3.3  Fields 
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Now  using  the  concept  of  the  group,  the  algebraic  structure  of  a  field  can  be  intro¬ 
duced.  A  field  ( F )  is  defined  such  that  it  can  be  thought  of  as  a  combination  of  a  com¬ 
mutative  group  under  addition  ( +  )  and  a  commutative  group  under  multiplication  ( ■  )• 
Basically,  this  makes  the  field  a  set  of  elements  in  which  the  operations  of  addition  (  + ), 
multiplication  (  •  ),  subtraction  ( - ),  and  division  (  4- )  are  defined.  The  follow  properties 
must  be  satisfied  to  define  a  field: 

a)  F  is  a  commutative  group  under  ( +  )  with  additive  identity  element  (  0 ). 

b)  The  set  of  nonzero  elements  in  F  is  a  commutative  group  under  ( * ) 
with  multiplicative  identity  element  ( 1 ). 

c)  Multiplication  ( • )  distributes  over  addition  ( + ),  that  is  V  a,b,c  e  F; 
a  •  (b  -I-  c)  =  (a  •  b)  +  (a  '  c) 

In  block  coding,  fields  of  finite  order  are  key  algebraic  structures  of  interest.  A  finite 
field  is  referred  to  as  a  Galois  field  that  is  named  after  Evariste  Galois,  the  discoverer  of 
finite  fields.  A  Galois  field  is  denoted  by  GF(q),  where  the  order  of  the  field  is  q. 

It  should  be  obvious  fi'om  the  definition  that  a  field  needs  to  contain  at  least  two  ele¬ 
ments,  the  additive  identity  element  ( 0 )  and  the  multiplicative  identity  element  ( 1 ). 

The  field  that  is  formed  from  these  two  elements  is  known  as  the  binary  field  GF(2).  The 
binary  field  and  its  extensions  are  very  important  and  have  many  applications  in  engi¬ 
neering.  Before  presenting  the  binary  extension  field,  primitive  elements  and  primitive 
polynomials  must  be  introduced.  But  first  the  concept  of  a  vector  space  will  be  intro¬ 
duced. 
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3.4  Vector  Spaces 


Let  V  represent  a  set  of  elements  called  vectors  and  F  represent  a  field  of  elements 
called  scalars.  Two  additional  operations  are  introduced  that  can  be  applied  to  the  ele¬ 
ments  of  V  and  F.  Let  (  +  )  be  a  binary  additive  operation  that  maps  linear  combinations 
of  vectors  v„  Vj,  Vj, ...,  contained  in  V  into  a  single  vector  V;  also  contained  in  V.  This 
binary  additive  operation  is  referred  to  as  vector  addition.  And  let  ( *  )  represent  a  binary 
multiplicative  operation  performed  on  a  scalar  aj  contained  in  F  and  a  vector  V;  contained 
in  V,  such  that  (  •  )  maps  a^  and  V;  into  a  vector  Uj  =  a;  •  V;  also  contain  in  V.  A  vector 
space  is  formed  if  the  following  rules  are  meet: 

a)  V  forms  a  commutative  group  under  ( + ). 

b)  V  a  e  F  and  Vv€V;a-v  =  u€V. 

c)  Associativity  under  multiplication  ( •  ):  V  a,b  e  F  and  V  v  e  V; 

(a  •  b)  •  V  =  a  •  (b  •  v) 

d)  Multiplication  (  •  )  distributes  over  addition  (  +  ):  V  a,b  e  F  and  V  u,v  e  F; 
a  •  (u  +  v)  =  (a  *  u)  +  (a  •  v),  and  (a  +  b)  •  u  =  (a  *  u)  +  (b  *  u). 

e)  V  V  €  V,  1  •  V  =  V  •  1  =  v;  the  scalar  multiplication  identity  element  is  1. 

In  a  vector  space,  F  is  commonly  referred  to  as  the  “ground  field”  or  “base  field”  of 
the  vector  space.  The  properties  of  the  vector  space  allow  for  the  existence  of  what  is 
known  as  a  vector  n-tuple.  A  n-tuple  has  the  form  v  =  (Vo,  v„  V2,  ...,v„.,)  where  the  ele¬ 
ments  (Vj)  are  from  the  ground  field  F.  It  is  the  n-tuple  representation  of  a  vector  space 
that  is  so  important  in  implementing  error  control  codes  due  to  the  association  of  binary 
words  with  vector  n-tuples.  For  example,  if  the  ground  field  is  taken  to  be  binary 
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( GF(2) ),  then  00101  is  a  five  bit  binary  v^ord  can  be  represented  by  v  =  (00101)  a  5- 
tuple  taken  from  the  vector  space  V5 .  Of  course,  here  the  vector  space  is  finite  and  com¬ 
posed  of  32  (2®)  vector  5-tuples. 

Some  other  important  ideas  of  vector  spaces  associated  with  error  control  coding  are 
the  dimension  of  the  vector  space,  the  basis  of  a  vector  space,  and  the  spaiming  set  of  a 
vector  space.  A  spanning  set  S  =  (vq,  v„  Vj,  ...,v„.,),  is  a  set  of  vectors  fi'om  which  the  lin¬ 
ear  combination  of  the  vectors  (ao’Vo  +  a,-v,  +  a2’V2  +...+  an.,'V„.i  ;  where  a^  =  0  or  1) 
forms  all  other  possible  vectors  of  the  vector  space.  Going  another  step  further,  a  basis  is 
a  set  of  vectors  of  minimum  order  that  spans  the  vector  space.  Therefore,  a  basis  is  a 
spanning  set  that  has  no  additional  redundant  vectors.  A  basis  set  is  always  a  spanning 
set  but  a  spanning  set  may  not  be  a  basis.  Finally,  the  dimension  of  a  vector  space  is 
equal  to  the  number  (k)  of  vectors  in  the  basis  set.  This  of  course  leads  to  the  fact  that 
each  vector  of  a  vector  space  is  xmique. 

3.5  Primitive  Element 

Elements  of  a  finite  field  have  individual  cardinalities.  Defining  a^  —  a  •  a,  a^  =  a  •  a  • 
a,  etc...,  the  order  of  an  element  (  n  =  ord(a) )  is  the  value  that  satisfies  a"  =  1 ,  the  identity 
element.  Of  particular  importance,  are  the  elements  in  field  GF(q)  with  order  q-1,  which 
are  called  primitive  elements.  It  can  be  proven  that,  since  a  primitive  element  has  order 
equal  to  one  less  than  that  of  the  size  of  the  field,  consecutive  powers  of  a  primitive  ele¬ 
ment  can  represent  all  nonzero  elements  of  a  field.  It  also  has  been  proven  that  there  exist 
at  least  one  element  in  GF(q)  with  order  q-1,  and  hence  there  is  always  at  least  one  primi- 
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tive  element  associated  with  a  finite  field.  The  importance  of  primitive  elements  will  be 
further  illustrated  in  the  section  on  the  development  of  binary  extension  fields. 

3.6  Primitive  Polynomials 

It  has  been  shown  thus  far  that  a  field  GF(q)  can  be  represented  by  0  and  q-1  consecu¬ 
tive  powers  of  a  primitive  element  (  a  ).  Remembering  the  additive  and  multiplicative 
nature  of  the  field,  a  field  can  be  composed  of  polynomials  of  m*  degree  whose  coeffi¬ 
cients  are  taken  from  a  defined  finite  field.  The  notation  GF(q)[x]  represents  such  a 
polynomial  f(x)  with  coefficients  taken  from  GF(q).  As  defined  earlier,  this  collection  of 
polynomials  is  called  a  commutative  ring  with  identity. 

A  polynomial  f(x)  which  cannot  be  factored  into  the  product  of  lower  degree  polyno¬ 
mials  in  GF(q)[x],  is  an  irreducible  polynomial  in  GF(q).  That  is  to  say,  the  polynomial 
f(x)  is  in  the  simplest  form  over  GF(q)[x].  For  example,  f(x)  =  x^  +  x  +  1  is  irreducible  in 
GF(2)[x],  however  f(x)  can  be  factored  in  GF(4)[x].  With  this  illustration,  it  is  clear  that 
when  defining  irreducibility  the  field  must  be  specified. 

An  irreducible  polynomial  p(x)  of  arbitrary  degree  m,  defined  in  GF(p)[x]  where  p  is  a 
prime,  is  called  a  primitive  polynomial  if  the  smallest  value  of  n  for  which  p(x)  divides 
x"  -  1  is  n  =  p"*  - 1 .  For  example,  p(x)  =  x^  +  x  +  1  is  a  primitive  polynomial  in  GF(2)[x] 
and  the  smallest  degree  polynomial  of  the  form  x"  - 1  that  p(x)  divides  is  x’  -  1 .  Note 
that  all  primitive  polynomials  in  GF(q)[x]  are  irreducible,  but  not  all  irreducible  polyno¬ 
mials  are  primitive.  An  extremely  important  property  of  primitive  polynomials,  which 
will  be  exploited  in  the  development  of  RS  codes,  is  that  the  roots  (a')  of  a  primitive 
polynomial  p(x)  in  GF(p)[x]  have  order  p""  -  1  .  Recall,  an  element  of  order  p™  -  1 
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forms  all  nonzero  elements  of  a  field  GF(q)  of  order  q  =  p"  - 1 .  Summarizing,  an  m*  de¬ 
gree  polynomial  in  GF(p)[x],  is  a  function  of  x  with  coefficients  in  GF(p)  and  has  roots 
(a')  over  the  field  GF(p"’ ). 

3.7  Binary  Extension  Fields 

In  the  study  of  error  control  coding  we  are  interested  in  polynomials  in  GF(p)[x]  that 
have  coefficients  in  GF(p)  and  roots  (a*)  in  GF(p‘”),  where  p  is  equal  to  the  binary  ground 
field.  Simply  stated,  p  =  2.  By  defining  a  field  in  this  fashion,  the  elements  of  the  field 
become  binary  n-tuples  of  length  m  over  GF(p'”).  Essentially,  binary  words  of  length  m 
are  defined.  Using  a  primitive  polynomial  p(x)  of  degree  m,  and  the  fact  a  is  a  primitive 
root  of  p(x),  the  construction  of  a  binary  extension  field  of  order  p”  will  now  be  illus¬ 
trated. 

a)  Step  1 :  Choose  a  primitive  polynomial  of  degree  m,  where  m  is  equal  to  the 
size  of  the  binary  word  of  interest.  For  instance,  if  working  with  3 -bit  binary 
symbols,  choose  a  primitive  polynomial  of  degree  three.  Choose  p(x)  =  x^  +  x 
+1. 

b)  Step  2:  Let  a  be  a  primitive  root  of  p(x).  Therefore,  p(a)  =  a^  +  a+l  =  0  and 

a^=  a  +  1. 

c)  Step  3:  Recall  that  p(x)  must  divide  f(x)  =  x^  -  1;  where  n  =  2^  -  1.  If  p(x)  di¬ 
vides  f(x)  and  has  a  as  a  root,  then  a  must  also  be  a  root  of  f(x).  Therefore, 

f(a)  =  a’  -  1  =  0  and  hence  a’  =  1 . 
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d)  Step  4:  Include  the  zero  element.  Since  the  consecutive  powers  of  the  primi¬ 
tive  element  a  only  form  the  nonzero  elements  of  the  field,  the  zero  element 
must  be  included. 

e)  Step  5:  Begin  construction  of  the  binary  extension  field  using  the  properties 
of  the  above  steps  to  keep  all  elements  in  their  reduced  format.  The  field  ele¬ 
ments  are 

0,  a°=l,  a'=a,  a^  =  a-i-l,  a'*  =  a^a  =  (a+1)  a  =  +  a, 

a*  =  =  (a^  +  a)a  =  (a^  +  a^)  =  (a  +  1  +  a^)  =  -i-  a  +  1,  etc... 


Table  3.0  Construction  of  Binary  Extension  Field  GF(8) 


EXPONENTIAL 

REPRESENTATION 

POLYNOMIAL 

REPRESENTATION 

BINARY  N-TUPLE 
REPRESENTATION 

0 

0 

000 

1 

1 

100 

a 

a 

010 

001 

a  +  1 

110 

a4 

+  a 

on 

+  a  H-  1 

111 

-f  1 

101 
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4  Reed-Solomon  Codes 


With  the  evolution  of  today’s  complex  communication  systems,  it  is  apparent  that  a 
strong  error-correcting  scheme  is  needed  to  combat  burst  errors.  Reed-Solomon  codes, 
once  thought  of  as  only  a  great  non-practical  mathematical  theory,  have  proven  to  be  the 
answer  to  the  burst  error-correcting  problem.  Reed-Solomon  codes  received  their  glory 
from  such  applications  as  the  deep  space  exploration  missions  of  Galileo  and  Voyager 
[McE2],  the  compact  disc  optical  recording  system  [Imm2],  and  many  applications  in 
magnetic  recording  systems  [Her].  With  the  advances  being  made  in  communication 
systems,  RS  codes  are  sure  to  see  many  future  applications.  The  remainder  of  this  section 
will  attempt  to  explain  what  RS  codes  are  and  how  they  are  encoded  and  decoded. 

Reed-Solomon  codes  are  a  subclass  of  linear  block  codes.  Like  any  other  block  codes, 
the  RS  codes  append  parity  information  to  a  block  of  data  so  that  at  the  receiver  error  de¬ 
tection  and  correction  is  possible.  Reed-Solomon  codes  are  typically  specified  by  their 
length  (n)  and  dimension  (k).  Therefore,  we  specify  a  particular  RS  code  as  a  (n,k)  RS 
code.  Recall  from  Section  3  that  k  is  the  number  of  elements  in  the  basis  of  the  vector 
space  and  is  an  indication  of  the  number  of  data  symbols  in  the  original  information 
block.  If  n  (where  n  >  k)  is  equal  to  the  number  of  symbols  in  the  encoded  block,  then 
there  is  n-k  parity  symbols  appended  to  the  k  original  data  symbols.  Also,  a  RS  code  is 
said  to  have  a  code  rate  defined  as  R^  =  k  /  n.  Refer  to  Figure  4.0  below  describing  the 
structure  of  a  (7,3)  RS  codeword. 
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n  symbols 


-► 


4 — m  bits — ^ 


n  -  k  parity  symbols - ►  •* - k  data  symbols 


Figure  4.0  Structure  of  a  (7,3)  Reed-Solomon  Code 


4.1  Error-correcting  Capability 

Reed-Solomon  codes  eire  word  oriented  linear  block  codes  capable  of  correcting  t  ran¬ 
dom  symbol  errors.  Reed-Solomon  codes  have  the  ability  to  correct  burst  errors  of  length 
equal  to  t  consecutive  symbols.  Since  RS  codes  are  word  oriented,  that  is  they  are  com¬ 
posed  of  m-bit  symbols,  they  can  correct  bursts  of  maximum  length  equal  to  t  •  m  bits.  It 
can  be  shown  that  the  guaranteed  length  for  burst  error  correction  of  a  RS  code  is 
b  =  m(t  - 1)  +  1  bits.  This  can  be  expanded  to  guaranteed  double  burst  error  correction 
(b  =  m(Lt/2j  -  1)  +  1)  and  guaranteed  triple  burst  error  correction  (b  =  m(Lt/3j  - 1)  +  1) 
[Kam2].  The  definitions  given  here  for  the  guaranteed  (not  maximum)  burst  error  cor¬ 
rection  will  be  utilized  in  the  following  sections  to  define  the  efficiencies  of  RS  codes. 


4.2  Maximum  Distance  Codes 

Reed-Solomon  codes  have  the  important  property  of  being  Maximum  Distance  Sepa¬ 
rable  (MDS)  codes.  It  has  been  shown  that  RS  codes  satisfy  the  Singleton  bound  with 
equality  and  therefore  are  considered  Maximum  Distance  codes  [Wicl].  Maximum  dis- 
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tance  codes  have  the  highest  obtainable  error  correction  capability  of  any  block  code  of 
the  same  redundancy  [Sin], 

Hamming  weight  has  been  defined  as  the  number  of  non-zero  elements  in  a  binary 
symbol.  For  example,  the  symbol  (1 1010001)  has  a  Hamming  weight  of  four  (denoted  W 
=  4).  Hamming  distance  is  defined  as  the  number  of  positions  between  two  binary  sym¬ 
bols  that  are  in  difference.  For  example,  the  Hamming  distance  between  the  two  binary 
symbols  (1 1010001)  and  (10001010)  is  equal  to  five  (denoted  D  =  5).  The  Hamming 
distance  is  also  shown  to  be  the  Hamming  weight  of  the  modulo-two  addition  of  any 
combination  of  symbols.  F or  example,  (1101 0001)  +  (1000 1010)  =  (01011011)  and  the 
Hamming  weight  of  the  resultant  (0101 101 1)  is  five  as  before.  The  minimum  distance 
(d^in)  of  a  code  is  defined  as  the  minimum  Hamming  distance  between  all  possible  valid 
pairs  of  codewords  in  a  given  vector  space.  Recall  that  a  vector  space  defined  over  a  bi¬ 
nary  extension  field  is  composed  2"  possible  code-vectors,  where  n  equals  the  code 
length.  For  the  same  vector  space,  there  are  2"^  valid  codewords,  where  k  is  the  dimension 
of  the  vector  space.  Therefore,  there  are  at  least  d„i„  bits  separating  valid  codewords. 
Maximum  distance  codes,  are  codes  that  have  the  largest  possible  minimum  distance  of 
separation  between  valid  codewords.  In  fact  it  has  been  shown  by  [Sin]  that  if  a  code  has 
d^i„  =  n  -  k  +  1 ,  it  is  a  MDS  code.  This  can  be  explained  by  visualization  of  any  two 
codewords  n  symbols  in  length  with  k  data  symbols  and  n-k  parity  symbols.  To  be  valid 
codewords,  the  k  data  symbols  must  at  minimum  differ  in  one  position.  Also,  the  most 
desired  case  would  be  when  all  n-k  parity  symbols  differ  between  two  valid  codewords. 
Combining  the  above  two  statements,  the  largest  minimum  separation  between  valid 
codewords  is  equal  to  the  maximum  possible  distance  between  all  parity  symbols 
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(D  =  n-k)  and  the  single  distance  between  data  symbols  ( D  =  1),  yields  d^jj^  =  n  -  k  +  1 . 


4.3  Encoding 

Since  the  introduction  of  RS  codes  in  1960  by  Irving  Reed  and  Gus  Solomon  [Ree], 
there  have  been  many  different  approaches  for  encoding  of  RS  codes.  These  differences 
are  due  to  the  fact  that  there  existed  no  efficient  methods  of  decoding  RS  codes,  therefore 
new  encoding  methods  where  created  that  gave  way  to  better  decoding  algorithms.  The 
original  method  defining  RS  codes  as  polynomials  over  certain  finite  fields  can  be  re¬ 
viewed  in  Reed  and  Solomon’s  original  1960  paper  [Ree].  There  is  also  a  summary  in 
[Wic2].  The  method  of  encoding  RS  codes  that  will  follow  is  the  extended  non-binary 
BCH  code  approach.  There  is  a  reason  for  choosing  the  extended  BCH  approach.  In 
most  applications  this  approach  leads  to  more  efficient  decoding  algorithms. 

In  order  to  construct  an  RS  code  as  an  extension  of  a  BCH  code,  one  must  use  the 
generator  polynomial  approach.  RS  codes  are  represented  by  generator  polynomials  g(x). 
These  polynomials  are  constructed  over  GF(2'")  and  therefore  have  as  roots  the  primitive 
elements  a' .  In  order  to  be  MDS,  g(x)  must  meet  the  BCH  bound.  Therefore,  the  a' 
primitive  elements  must  be  2t  =  n-k  consecutive  primitive  roots  of  g(x).  Such  a  polyno¬ 
mial  must  be  of  the  form  g(x)  =  (x  -  a'’)  (x  -  a*’*' )  (x  -  )...  (x  -  ).  For  narrow- 

sense  RS  codes,  g(x)  can  be  represented  by  the  following  equation  (in  which  b  equals 
one): 

2t 

g(x)=n  (x-a‘) 

i=l 
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When  all  the  factors  of  g(x)  are  multiplied  out  to  form  the  generator  polynomial,  it  will  be 
a  polynomial  of  degree  2t  =  n  -  k  defined  over  the  field  GF(q)  where  q  =  2"’  and  will  thus 
have  q‘‘  valid  codewords. 

By  defining  RS  codes  as  extended  BCH  codes  using  the  polynomial  approach,  the  cy¬ 
clic  properties  can  be  utilized.  The  generator  polynomial  has  the  form 
g(x)  =  go  +  giX  +  g2X^  + ...  +  gn-k  x”'‘‘ ,  and  a  block  of  information  data  can  be  represented 
as  a  polynomial  of  the  form  d(x)  =  do  +  d,x  +  d2x'  + ...  +  d^.,  x"'' .  Therefore,  a  codeword 
C  is  encoded  as  c(x)  =  d(x)g(x)  and  is  represented  as  c(x)  =  Cq  +  c,x  +  C2X  + ...  +  c„_jX  . 
Typically  RS  codes  are  put  into  a  systematic  format  in  order  to  reduce  the  number  of 
steps  needed  for  encoding.  A  codeword  is  in  systematic  format  if  it  has  the  form 
C  =  (Po,  Pi,  P2,  -  P2t.i ,  do,  d„  d2, ...  dk., )  ,  where  the  data  symbols  (d;)  are  grouped  to¬ 
gether  and  the  n-k  parity  symbols  (pj)  are  also  grouped  together.  Using  the  cyclic  proper¬ 
ties  and  polynomial  representations,  the  codeword  has  the  form  c(x)  =  p(x)  +  x  d(x)  , 
where  multiplying  by  x"''^  simply  performs  a  cyclic  shift  of  the  data  symbols  into  their 
proper  position.  In  the  above  polynomial  representation  the  parity  information  is 
equivalent  to  the  remainder  of  x"  M(x)/g(x),  that  is  p(x)  =  x-'MCx)  mod  g(x)  [Has].  The 
resultant  code  word  pol5momial  c(x),  systematic  or  non-systematic,  has  as  roots  the  same 
2t  consecutive  powers  of  cc  as  does  the  generator  polynomial.  This  property  allows  for  an 
easy  way  of  determining  if  a  valid  codeword  has  been  received  and  for  efficient  decoding 
algorithms  of  RS  codes. 
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4.4  Decoding 

Since  the  introduction  of  RS  codes  in  1960,  there  has  been  a  rush  by  researchers  to 
find  an  efficient  decoding  algorithm  for  these  codes.  The  simple  techniques  used  in  lin¬ 
ear  block  codes  just  do  not  work.  For  instance  using  the  syndrome  look-up  table  method 
is  incomprehensible  due  to  the  number  of  syndrome  computations  needed.  For  instance, 
a  (63,53)  RS  code  requires  approximately  10^°  syndrome  calculations  and  has  approxi¬ 
mately  10®^  (64”)  valid  codewords. 

In  the  original  paper  by  Reed  and  Solomon  [Ree]  they  introduced  a  possible  method 
for  decoding.  This  decoding  method  involved  the  solution  of  a  set  of  k  equations  in  k 
unknowns.  The  method  proved  to  only  be  useful  for  the  simplest  RS  codes.  Shortly  af¬ 
ter,  Peterson  [Pet2]  presented  a  popular  decoding  method  for  BCH  codes.  Gorenstein  and 
Zierler  [Gor],  and  later  Chien  [Chi]  and  Forney  [For2],  expanded  the  Peterson  algorithm 
for  BCH  codes  to  the  non-binary  case  for  RS  codes.  These  codes,  although  much  more 
efficient  than  any  others,  were  only  applicable  to  codes  with  a  small  number  of  errors. 

It  was  not  imtil  Berlekamp  [Ber2]  presented  his  version  of  an  efficient  decoder  for  use 
with  non-binary  BCH  and  RS  codes  that  the  application  of  RS  codes  became  more  than  a 
mathematical  exercise.  One  year  later  Massey  [Mas]  showed  that  the  Berlekamp  algo¬ 
rithm  could  be  realized  as  a  technique  of  finding  the  shortest  Linear  Feedback  Shift  Reg¬ 
ister  (LFSR)  capable  of  generating  a  given  sequence.  The  combination  of  the  two  ideas 
has  now  become  what  is  known  as  the  Berlekamp-Massey  algorithm.  Although  the  Ber- 
lekamp-Massey  algorithm  is  not  the  most  mathematically  simple  to  understand,  it  pro¬ 
duces  an  efficient  technique  for  decoding  RS  codes.  The  explanation  of  the  Berlekamp- 
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Massey  algorithm  that  follows  is  a  simple  overview  and  the  reader  is  referred  to  [Berl], 
[Ber3],  [Has],  [Lin],  [Nas],  [Wicl], 

The  Berlekamp-Massey  algorithm  is  a  four-step  process  for  error  correction  in  RS 
codes.  Reed-Solomon  codes,  which  unlike  BCH  codes  are  word  oriented,  require  an  ad¬ 
ditional  step  of  determining  the  error  value.  The  first  step  in  the  Berlekamp-Massey  algo¬ 
rithm  is  to  calculate  the  syndromes.  Next  the  error  locator  polynomial  must  be  obtained. 
Once  the  error  locator  polynomial  has  been  obtained,  the  error  location  can  be  found.  Fi¬ 
nally,  the  error  values  are  calculated.  These  four  steps  are  explained  in  the  following 
paragraphs. 

Recall  that  a  received  vector  has  the  form  r(x)  =  ro  +  rjX  +  r2X^  -i- ...  +  r„.,x"‘'.  The  re¬ 
ceived  vector  can  be  represented  has  the  combination  of  a  valid  codeword  and  an  error 
pattern,  r(x)  =  c(x)  -i-  e(x)  where  c(x)  =  Cq  +  c,x  +  CjX^  +  ...  +  c„_]X"  '  is  a  valid  codeword 
and  e(x)  =  Oq  -i-  e,x  +  ejX^  +  ...  +  e„.,x"''  is  an  error  pattern.  Since  the  codeword  c(x)  has  as 
roots  the  2t  consecutive  powers  of  a,  the  received  vector  will  have  the  same  roots  if  no 
errors  are  present.  That  is  r(x)  =  c(x)  when  e(x)  is  equal  to  all  zeros.  Calculating  the  syn¬ 
dromes  is  done  by  evaluating  the  received  vector  at  all  2t  consecutive  roots  of  the  code¬ 
words,  Sj  =  r(a‘)  =  c(a‘)  +  e(a‘)  for  i=  l,2,...,2t.  Since  S;  will  equal  zero  except  where 
e(a')  is  not  zero,  the  syndrome  is  calculated  from  the  equation  below  to  find  out  if  errors 


exist. 


The  second  step  of  the  Berlekamp-Massey  algorithm  is  to  obtain  the  error  locator 
polynomial  A(x).  As  the  name  implies,  this  polynomial  gives  information  as  to  the  loca¬ 
tions  of  the  errors  in  the  received  code- vector.  It  is  actually  the  integer  power  of  the  in¬ 
verse  root  of  the  locator  polynomial  that  reveals  the  error  location.  Simply  stated:  find 
the  roots  of  the  error  locator  polynomial,  find  the  inverse  of  the  roots,  and  then  the  power 
of  the  resultant  is  the  location  of  the  error.  For  example  if  the  root  was  o?  defined  over 
GF(8),  then  the  inverse  of  is  and  the  error  is  at  the  fifth  location. 

It  is  the  method  of  finding  the  error  locator  polynomial  that  makes  the  Berlekamp- 
Massey  algorithm  so  complex.  The  error  locator  polynomial  is  foimd  through  an  iterative 
process,  usually  performed  in  tabular  form,  whereby  the  solution  for  A(x)  is  checked  for 
correctness  after  each  iteration.  If  the  error  locator  polynomial  A(x)  is  not  correct,  then  a 
correction  factor  is  added  to  A(x)  and  the  process  continues  for  2t  iterations.  The  error 
locator  polynomial  A(x)  is  given  below  for  v  errors,  where  Aj  is  a  coefficient  of  the  form 
aT 

A(x)  =  (l+Ai  X+Aj  X^+...+Av  X'') 


After  the  error  locator  polynomial  A(x)  has  been  found,  the  third  step  is  to  find  all  of 
the  roots  of  the  locator  polynomial.  This  can  be  performed  by  factoring  A(x)  into  the 
form  A(x)  =  (l+AjX’Xl+Ai+iX'^')...  (1+A^,X')  or  just  by  evaluating  A(x)  at  all  power  of 
a.  Once  the  roots  of  A(x)  are  found,  the  inverses  of  the  roots  reveal  the  error  locations. 

The  fourth  step  in  the  Berlekamp-Massey  algorithm  is  that  of  calculating  the  actual 
error  values  at  the  now  known  error  locations.  In  order  to  perform  this  task,  the  polyno- 
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mial  known  as  the  error  magnitude  polynomial  Z(x)  is  defined  using  the  error  locator 
polynomial  coefficients  and  the  syndrome  values.  The  error  magnitude  polynomial  is 
defined  as  Z(x)  =  1  +  (S,  +  A,)x  +  (S,  +  A,  S,  +  A2)  x'  +...+  (S,  +  A,  S,.,  +...+  A,)  x"  for 
V  errors.  This  error  magnitude  polynomial  is  used  in  the  following  equation  to  find  the 
error  values.  The  munerator  of  the  error  value  equation  is  the  error 

Z(X-‘) 

e..  = - ^ - - - 

n  (i+XiX-^) 

magnitude  polynomial  evaluated  at  the  root  of  the  error  locator  polynomial.  The  de¬ 
nominator  is  equal  to  v  factors  where  Xj  represents  all  roots  of  the  error  locator  polyno¬ 
mial  that  are  not  at  the  current  location  being  calculated,  and  X' j  is  the  inverse  root  of  the 
error  locator  polynomial  for  the  current  error  value.  Once  all  v  values  for  Oy  are  foimd  the 
error  polynomial  can  be  constructed  and  added  back  to  the  received  vector  for  error  cor¬ 
rection.  The  error  polynomial  would  have  the  form  e(x)  =  q  +  ej+ix"^'  +  ei+2x'^^  +  ...  +  e„. 
,x"''  and  the  corrected  codeword  are  obtained  by  c(x)  =  r(x)  +  e(x). 

5  Interleaving 

Communication  systems  discussed  up  to  this  point  have  been  shown  to  operate  over 
channels  with  memory.  Furthermore,  it  is  the  memory  in  the  channel  that  is  to  blame  for 
the  majority  of  burst  errors.  Interleaving  is  a  technique  that  can  be  applied  in  an  attempt 
to  randomize  the  transmission  of  encoded  signals  over  the  bursty  chaimel.  This  random¬ 
izing  of  the  data  allows  the  channel  to  appear  memoryless.  The  remainder  of  this  section 
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will  discuss  the  two  most  widely  applied  interleaving  techniques,  block  interleaving  and 
convolutional  interleaving. 

5.1  Block  Interleaving 

Block  interleaving  attempts  to  scramble  the  codeword  over  the  chaimel  so  that  the 
channel  appears  to  be  random.  Block  interleaving,  although  appearing  random  to  the 
channel,  is  truly  a  one-to-one  deterministic  ordering  of  the  codeword  symbols  [Ram]. 

This  one-to-one  sequencing  of  the  data  provides  for  a  very  efficient  process  in  that  un¬ 
scrambling  of  codewords  is  easily  earned  out.  The  unscrambling,  or  deinterleaving  as  the 
process  will  be  called  hereafter,  is  essentially  an  interleaver  used  in  reverse  order.  All  of 
the  important  properties  will  be  brought  forth  in  the  following  discussion  of  the  block  in¬ 
terleaving  process. 

A  block  interleaver  is  nothing  more  than  a  buffering  mechanism  inserted  between  the 
error  coding  circuit  (channel  encoder)  and  the  modulator.  Similarly,  the  de-interleaver  is 
a  buffer  placed  after  the  demodulator  and  before  the  channel  decoder.  Figure  5.0  is  a 
simplified  block  diagram  of  a  typical  communication  system.  The  interleaver  and  de¬ 
interleaver  can  be  thought  of  as  a  d  x  n  matrix  where  d  is  the  depth  of  the  interleaver  indi¬ 
cating  the  number  of  codewords  stored.  Each  cell  in  the  matrix  represents  one  s5unbol  of 
length  m  bits.  Depending  on  the  size  of  the  interleaver  desired,  it  is  normally  imple¬ 
mented  with  memory  devices.  Therefore,  in  the  case  of  RS  codes,  memory  devices  of 
n-m-d  bits  are  required. 
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Figure  5.0  Simplified  Communications  System  Block  Diagram 

The  operation  of  the  block  interleaving  process,  as  illustrated  m  Figure  5.1,  is  simply 
the  inputting  of  codewords  into  the  interleaver  matrix  in  a  row-by-row  fashion,  then  out- 
putting  them  to  the  modulator  column-by-column.  This  technique  acts  to  spread  the 
codewords  over  the  channel  and  therefore  spreads  the  effects  of  a  burst  error  over  multi¬ 
ple  codewords.  With  block  interleaving,  as  shown  in  Figure  5.2,  a  burst  error  of  length  b 
<  d  in  the  binary  case  or  b  <  m  d  in  the  non-binary  case,  will  be  reduced  to  a  single  er¬ 
ror  in  each  codeword  [Ste].  Therefore,  the  length  of  the  burst  error  to  be  corrected  is  de¬ 
pendent  on  the  depth  (d)  of  the  block  interleaver.  The  depth  of  the  block  interleaver  can 
be  increased  but  not  without  the  tradeoff  of  increased  delays.  Also  as  shown  in  Figure 
5.3,  any  burst  error  of  length  b  =  r-m-d  where  r  >  1,  will  produce  r  symbol  errors  in  each 
codeword.  Hence,  a  code  with  error-correcting  ability  t  >  r  would  be  needed.  Therefore, 
by  using  block  interleaving  techniques,  a  less  complex  of  the  error  correction  code  can  be 
applied.  Reed-Solomon  codes  with  lower  error-correcting  capabilities  can  be  used  in 
conjunction  with  interleaving  techniques  to  correct  large  burst  errors. 
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Output  Stream:  A1 31 ,  Cl ,  D1 ,  El ,  FI .  A2,  B2,  C2,  D2,  E2,  F2,  A3, B3,  C3.  D3,  E3.  F3,  A4 .  C12,  D12,  E12,  FI2 

Figure  5.1  Block  Interleaver 

Conversely,  using  a  RS  code  with  a  higher  error  correction  capability  and  applying  inter¬ 
leaving,  further  strengthens  the  code  and  allows  for  the  correction  of  much  longer  burst. 
Note,  all  data  in  the  following  figures  are  read  from  left  to  right  with  the  left  most  data 
element  being  the  first  element  input  to  the  interleaver/de-interleaver.  The  shaded  rectan¬ 
gles  represent  the  occurrence  of  symbol  errors. 
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Input  Stream:  Al,  A2,  A3,  A4,  A5,  A6,  A7,  A8,  A9,  AIO,  A1 1,  A1 2,  B1 ,  B2,  B3,  B4,  B5,  B6, . ,F9,F10,  FlI ,  F12 


n  symbols  - ^ 


Al 

"A2 

A4 

~ 

A6 

~ 

A8 

A9 

AlO 

All 

A12 

■TT 

"bT 

1b7 

“b7 

Bll 

BI2 

Cl 

TT 

C4 

TT 

C6 

"c7 

C8 

C9 

CIO 

Cll 

C12 

TT 

“57 

"ST 

"57 

mo 

dTT 

577 

TT 

“ 

TT 

E6 

"e7 

"eT 

"e7 

eTT7 

ei7 

eTT 

FI 

"fT 

TT 

T? 

TT 

TT 

Fk 

"f7 

fT^ 

Fll 

fiT 

Output  Stream;  Al,Bl,Cl,DlJ 


E2,  F2,  A3,  B3,  C3,  D3,  E3,  F3,  A4, . ,C12,  D12,  E12,  F12 


burst 


De-inter  leaver 


Input  Stream:  A1,B1,C1,D1, 


E2,  F2,  A3,  B3,  C3,  D3,  E3,  F3,  A4, . ,C12,  D12,  E12,  F12 


n  symbols 


Al 

"aT 

A4 

A5 

A6 

"aT 

A8 

A9 

AlO 

All 

A12 

"bT 

w 

"bT 

■bT 

■s? 

“ 

isr 

B9 

bT7 

Bll 

B12 

"cT 

E 

“ 

"cT" 

C5 

C6 

"C7" 

C8 

C9 

CIO 

Cll 

C12 

"ST 

“ 

d4 

"57 

“57 

"57 

515 

5TT 

D12 

w 

E2 

"eT 

"eT 

"eT 

“ 

"eT 

"e7 

eTo 

eTT 

eiT 

W 

?"  .Jv 

“T7 

"fT 

"FT 

TT 

F6 

TT 

“TT 

"f7 

FiTT 

Fll 

fTT 

Output  Stream:  A1 


A3,  A4,  A5,  A6,  A7,  A8,  A9,  AI 0,  A1 1 ,  A1 2,  B 1 


m 


B3,  B4,  B5,  B6, . ,F9,F10,  FI  1,  F12 


Figure  5.2  Block  Interleaver  with  (b  <  d) 
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Figure  5.3  Block  Interleaver  with  (b  =  rd;  r  =  2) 
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5.2  Convolutional  Interleaving 

Convolutional  interleavers  are  very  much  like  block  interleavers  in  that  their  function 
is  to  scramble  the  codewords  before  transmission.  The  scrambling  of  a  convolutional  in¬ 
terleaver  is  also  a  one-to-one  deterministic  process.  It  can  be  seen  that  for  a  block  inter¬ 
leaver,  after  a  column  of  d  code  symbols  has  been  output  from  the  interleaver,  storage 
space  is  still  allocated  for  the  already  transmitted  column  of  data.  The  additional  storage 
needed  for  the  block  interleaver  produces  unwanted  delays  and  memory  requirements  in 
the  communication  system.  Convolutional  interleaving  is  a  more  optimal  interleaving 
technique  in  that  it  reduces  the  amount  of  storage  space  required  and  therefore  reduces 
delays.  There  are  many  excellent  references  available  describing  interleaving  techniques, 
however  the  explanation  of  convolutional  interleaving  that  follows  is  as  proposed  by 
[Fori]  and  further  explained  by  [Ste]  and  [Siml]. 

Figure  5.4  is  a  representation  of  a  convolutional  interleaver  and  corresponding  de¬ 
interleaver.  The  first  input  line  of  the  interleaver  is  connected  directly  to  the  channel 
(modulation  is  assumed).  All  other  input  lines  are  composed  of  registers  with  each  suc¬ 
cessive  register  containing  i  additional  storage  cells.  A  serial  data  stream  of  codewords 
from  the  encoder  circuit  is  input  to  the  convolutional  interleaver.  Serial  to  parallel  con¬ 
version  maybe  needed.  The  second  input  line  is  input  to  a  shift  register  of  length  I  sym¬ 
bols.  The  third  line  is  input  into  a  shift  register  of  length  2  f  symbols  and  the  process 
continues  for  d-1  lines,  for  a  total  of  d  input  lines.  Remember  that  for  RS  codes  each 
symbol  is  m  bits  in  length  and  therefore,  each  shift  register  cell  is  considered  to  be  m  bits 
in  length.  At  the  same  time  that  the  input  demultiplexer  is  inputting  a  codeword  symbol. 


292 


the  output  multiplexer  is  selecting  the  oldest  symbol  shifted  out  from  each  shift  register. 
The  output  stream  is  then  sent  to  the  modulator. 


Figure  5.4  Convolutional  Interleaver 


For  the  case  of  using  convolutional  interleavers  to  scramble  RS  coded  information,  it  is 
convenient  to  design  the  interleaver  with  the  depth  (d)  equal  to  the  length  of  a  codeword 
(d  =  n).  Therefore,  a  delay  of  d  ^  symbols  or  d  bits  is  placed  between  adjacent  input 
symbols.  From  this  implementation,  it  is  obvious  that  a  burst  of  length  (d  ^+l)m  bits  will 
only  introduce  one  symbol  error  in  d  ^  consecutive  codewords  after  de-interleaving  takes 
place. 


5.3  Pseudorandom  Interleaving 

Another  important  interleaving  technique  worth  mentioning  is  called  pseudorandom 
interleaving.  Pseudorandom  interleavers  work  on  the  principle  of  using  synchronous 
pseudorandom  number  generators  in  both  the  interleaver  and  deinterleaver.  The  pseu¬ 
dorandom  codes  are  used  to  vary  the  depth  of  interleaving  in  a  manner  only  known  to  the 
communication  system.  The  importance  of  this  technique  is  in  combating  intentional 
jamming  in  hostile  communication  environments.  If  the  depth  of  an  interleaver  is  known, 
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a  pulsed  jammer  can  be  used  with  a  pulse  repetition  interval  (PRI)  equal  to  the  interleav¬ 
ing  depth  in  order  to  produce  burst  errors  at  the  output  of  the  deinterleaver  [Sim2].  By 
using  pseudorandom  interleaving  depths,  pulsed  jammers  lose  their  effectiveness.  The 
information  provided  here  on  pseudorandom  interleaving  is  only  for  completeness  and 
will  not  be  explored  any  further. 

5.4  Interleaving  Delay 

Referring  back  to  Figure  5.1,  the  block  interleaving  scheme  uses  d  *  n  symbol  storage 
cells  or  d  •  m  •  n  bit  storage  cells.  This  introduces  a  delay  of  n(d-l)+l  symbols  because 
transmission  can  begin  as  soon  as  the  first  symbol  in  the  last  row  is  input.  Obviously, 
this  allows  for  the  output  of  the  first  column  before  the  input  of  the  remaining  symbols  in 
the  last  row.  If  measuring  the  delay  of  the  interleaver  in  bits,  then  there  is  a  m[n(d-l)+l] 
bit  delay.  The  delay  can  be  measured  in  time  (seconds)  using  the  known  fact  that  the  data 
rate  (  R  )  is  equal  to  the  inverse  of  the  bit  transmission  rate  (R  -  l/tg).  Therefore,  the  time 
delay  is  (ts'Tg)  seconds  where  Tg  =  m[n(d-l)+l]  is  the  interleaver  bit  delay  measured  by 
the  number  of  bits.  Therefore,  the  end-to-end  delay  time  for  block  interleaving  is 
2(ts’TB)  +  channel  delays. 

Convolutional  interleavers  are  more  efficient  than  block  interleavers  because  of  their 
smaller  storage  requirement  and  therefore  smaller  delay  times.  Referring  to  Figure  5.4,  it 
can  be  seen  that  the  total  number  of  storage  cells  needed  for  a  convolutional  interleaver  is 

n*  f(d-l)/2  symbols.  For  simplicity  define  I  to  be  equal  to  one  and  recall  each  symbol  is 
m  bits  in  length.  Therefore,  the  storage  requirement  can  be  measured  as  m[n(d-l)]/2  bits. 
As  it  was  done  for  the  block  interleaving,  time  delay  can  be  measmed  as  ts*m*[n(d-l)]/2 
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seconds.  Recalling  that  for  the  block  interleaver  the  storage  requirement  is  n(d-l)+l 
symbols  and  n(d-l)/2  symbols  for  the  convolutional  interleaver,  therefore  it  can  be  seen 
that  the  convolutional  interleaver  is  nearly  twice  as  efficient  as  the  block  interleaver.  The 
interleaver  efficiency  is  examined  in  more  detail  in  [Wicl]. 
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6  Interleaved  Burst  Error  Correction  Efficiency 

Up  to  this  point  a  preliminary  review  of  RS  codes  and  error  control  coding  in  general 

have  been  presented.  This  section,  along  with  Section  7,  constitutes  the  findings  of  the 
research  performed  on  evaluating  the  burst  error  correction  of  RS  codes.  First  in  this 
section  the  burst  error  correction  efficiency  of  linear  block  codes  will  be  explained  using 
the  traditional  definition  of  efficiency.  The  concept  of  phased  burst  error  correction  effi¬ 
ciency  will  also  be  explained.  Then  finally  the  burst  error  correction  efficiency  of  RS 
codes  will  be  defined. 

Recall  from  Section  2  one  that  burst  errors  are  defined  according  to  the  length  of  the 
burst  b,  and  the  guard  space  g  separating  the  individual  burst  errors.  Hence,  the  Gallager 
bound  states  that  for  a  binary  code  of  burst  error-correcting  capability  b  and  guard  space 
g,  the  code  rate  R  must  satisfy  the  inequality  below  [Gal2].  The  Reiger  bound,  which  is 


derived  from  the  Gallager  bound,  states  that  a  (n,k)  linear  block  code  with  burst  error- 
correcting  capability  b  has  to  meet  the  following  condition  [Cap]. 

n  —  k  >2b 

Also,  a  linear  block  code  with  burst  error  correction  b  and  detection  d  must  meet  the  con¬ 
dition  below. 

n  —  k  +  \>  b  +  d 
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Based  on  the  Reiger  bound,  the  traditional  method  for  defining  burst  error  correction  effi¬ 
ciency  is  z  as  shown  below.  This  ratio  is  called  the  “phased  burst  error  correction  effi¬ 
ciency”.  Linear  block  codes  for  which  z  =  1  are  considered  optimal  codes  because 


these  give  the  most  error  correction  achievable  for  the  amount  of  parity  used. 

6. 1  Maximum  Burst  Error  Correction  of  Reed-Solomon  Codes 

The  Gallager  boimd,  Reiger  Bound,  and  the  phased  burst  error  correction  efficiency 
are  defined  for  binary  linear  block  codes  correcting  single  burst  errors.  Reed-Solomon 
codes  are  non-binary  multiple  burst  error-correcting  linear  block  codes.  Recall  from  Sec¬ 
tion  4  that  RS  codes  have  random  error-correcting  capability  t  =  (n-k)/2  and  maximum 
single  burst  error-correcting  capability  b,max=  hn,  where  the  subscript  (1)  indicates  single 
burst  error  correction.  The  term  maximum  is  used  because  the  burst  errors  must  be  lo¬ 
cated  exactly  over  t  symbols  (which  may  not  always  be  the  case).  This  (bi^aJ  fonns  an 
upper  bound  on  the  single  burst  error-correcting  capability  of  RS  codes.  Using  the  phased 
single  burst  error-correcting  capability  of  RS  codes  b,n,ax  1  “  (n-k)/2,  it  can  be  seen 

that  RS  codes  are  always  optimal  phased  single  burst  error-correcting  codes. 
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6.2  Multiple  Burst  Error  Correction  of  Reed-Solomon  Codes 


In  some  situations  frequent  bursts  occur  in  a  given  codeword  and  single  burst  error- 
correcting  codes  are  of  no  use.  This  usually  occurs  when  the  guard  spaces  are  not  well 
defined  and  when  a  cluster  of  bursts  are  too  long  to  be  treated  as  a  single  large  burst.  As 
mentioned  earlier  RS  codes  are  multiple  burst  error-correcting  codes,  which  has  made 
them  very  popular  in  combating  these  types  of  errors.  With  this  in  mind,  the  phased 
(maximum)  double  and  triple  burst  error-correcting  capability  of  RS  codes  are  given  be¬ 
low.  Where  the  function  W  is  the  floor  function  and  is  defined  as  the  largest  integer  less 
than  or  equal  to  x.  Because  RS  codes  are  word-oriented  codes,  the  term  phased  burst 

.  1  M 

^  2  m  a  X 

^  3  m  ax 

error  correction  applies  only  to  the  case  that  the  burst  errors  occur  over  exactly  t  symbols. 
Therefore,  the  phased  burst  error  correction  of  RS  codes  is  the  maximum  burst  error  cor¬ 
rection. 

The  maximum  double-burst  error  correction  capability  is  defined  as  the  maximum  size 
of  two  simultaneous  burst  errors  that  can  be  corrected.  The  code  may  be  able  to  correct 
different  combinations  of  double  burst  errors  so  long  as  t  is  not  exceeded.  However,  no 
two  burst  greater  than  b2,„ax  corrected  at  once.  For  instance  if  m  =  8  and  t  =  6 

((255,243)  RS  code),  the  code  will  correct  two  burst  of  maximum  length  24  bits.  How¬ 
ever,  the  same  code  may  correct  two  bursts,  one  of  length  16  bits  and  the  other  32  bits.  If 
one  biurst  is  of  length  greater  than  b2n,ax5  then  the  other  burst  will  be  less  than  b2max-  The 
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maximum  burst-error  correction  capability  also  applies  to  triple  and  higher  order  burst 
error-correcting  codes.  The  general  form  for  phased  burst  error  correction  efficiency  of  a 
RS  codes  is  given  as  below.  It  can  be  easily  seen  that  although  the  double,  triple,  and 
higher  order  phased  burst  error  correction  efficiencies  of  RS  codes  are  not  optimal,  they 
have  very  high  efficiencies. 


where; 


2sb^ 
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6.3  Guaranteed  Burst  Error  Correction  Of  Reed-Solomon  Codes 

The  maximum  or  phased  burst  error  correction  capabilities  defined  above  are  only  ap¬ 
plicable  when  the  burst  errors  occur  exactly  over  t  symbols.  In  most  practical  applica¬ 
tions  the  burst  errors  will  not  occur  exactly  over  t  symbols  and  therefore  the  need  to  find 
the  guaranteed  burst  error-correcting  capabilities  of  a  code.  Recall  from  Section  4  that 
the  guaranteed  single,  double,  and  triple  burst  error  correction  capability  of  RS  codes  are 
b,,  \>2,  and  b3  restated  below.  This  is  the  burst  error-correcting  capability  of  a  RS  code  no 
matter  where  the  burst  errors  occur. 

6,  =  m  (t  -  1)  +  1 
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At  this  point  the  guaranteed  burst  error-correcting  efficiency  for  a  given  RS  code  is 
defined.  Using  the  traditional  definition  of  efficiency  and  the  guaranteed  burst  error- 
correcting  capability  equations,  the  single,  double,  and  triple  guaranteed  burst  error- 
correcting  efficiency  of  RS  codes  are  y„  Y2_andY3  as  follows. 

2[?»(?-l)  +  l] 

m(n-k) 


Ji  = 


r  ( 

t 

>1 

4 

.2. 

-1 

J 

+  1 

m(n  —  k) 


r  ( 

t 

>1 

6 

.3. 

-1 

J 

+  1 

m(n  -  k) 


For  the  general  case,  the  burst  error  correction  efficiency 
code  is  given  [KamS]. 


of  an  s  burst  error-correcting  RS 


where; 


m(n  —  k) 


/ 

t 

=  m 

\ 

y 

6.4  Analysis 

Reed-Solomon  codes  of  length  n  =  255,  symbol  length  m  =  8  bits,  are  very  popular 
and  find  many  applications  due  to  their  byte  sized  symbols.  Also,  codes  with  rate  R  less 
than  fifty  percent  are  not  often  used  because  of  the  large  amount  of  redundancy.  There¬ 
fore,  Table  6.0  shows  RS  single,  double,  and  triple  guaranteed  burst  error  correction  ca- 
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pabilities  and  efficiencies  for  codes  of  length  n  =  255  and  rate  R  >  50%.  For  a  complete 
table  of  RS  codes  of  length  n  =  7, 15, 31, 63, 127,  and  255  see  appendix  A. 

As  can  be  seen  jfrom  Table  6.0  and  as  illustrated  in  Figure  6.0,  the  guaranteed  single 
burst  error  correction  efficiency  approaches  optimality  very  quickly  as  code  rate  de¬ 
creases.  Note  that  as  can  be  observed  fi'om  the  equations  for  guaranteed  burst  error  cor¬ 
rection  efficiency  of  RS  codes,  Yi ,  Y2 ,  and  YsCan  never  equal  one,  however  they  are  re¬ 
markably  close.  The  fluctuations  in  Y2  and  Y3  are  due  to  the  fact  that  there  are  RS  codes 
that  have  the  same  double  and  triple  guaranteed  burst  error-correcting  capabilities  with 
increased  redimdancy  as  observed  in  the  flat  parts  in  the  graphs  of  Figure  6.1.  Figure  6.1 
only  shows  RS  codes  up  to  a  code  rate  of  eighty  percent  for  better  resolution.  For  a  com¬ 
plete  set  of  graphs  of  burst  error-correcting  capabilities  of  RS  codes  see  appendix  A. 

Using  the  guaranteed  burst  error-correcting  capability  and  efficiency  of  RS  codes,  a 
code  designer  now  has  additional  information  at  his/her  disposal  for  selecting  the  best 
possible  RS  code  for  the  error  correction  scheme  being  designed.  For  instance,  if  the  de¬ 
signer  is  interested  in  double  burst  error-correcting  codes,  the  codes  on  the  flat  areas  of 
the  guaranteed  double  burst  error-correcting  capability  graph  with  the  highest  code  rate 
are  the  best  choice.  The  other  codes  on  each  flat  region  have  decreased  code  rates  and 
lower  burst  error-correcting  efficiencies.  The  information  presented  here  can  help  guide 
the  code  designer  to  get  more  error  correction  for  the  amount  of  redundancy  (more  bang 
for  the  buck). 
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Table  6.0  Burst  Error  Correction  of  Reed-Solomon  Codes  of  Length  n  —  255 


(n.k)  t 
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Burst  Error  Correcting  Capability  (bits)  EFFICIENCY 


BURST  ERROR  CORRECTION  EFFICIENCY 
(Reed-Solomon  Codes  of  Length  n  =  255) 


—•—Single  Burst 
-ah- Double  Burst 
—•—Triple  Burst 


Figure  6.0  Guaranteed  Burst  Error  Correction  Efficiency 


GUARANTEED  BURST  ERROR  CORRECTING  CAPABILITY 
(Reed-Solomon  Codes  of  Length  n  =  255) 


Figure  6.1  Guaranteed  Burst  Error  Correction  Capability 
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7  Bxirst  Error  Correction  Efficiency 

Up  to  this  point  the  classical  definition  of  efficiency  for  a  linear  block  code  has  been 
expanded  to  include  non-binary  RS  codes.  This  research  has  produced  efficiency  equa¬ 
tions  that  will  be  useful  in  designing  RS  coding  schemes.  However,  due  to  the  fact  that 
interleaving  is  commonly  used  in  combating  bursty  errors,  a  look  at  the  efficiency  of  in¬ 
terleaved  RS  codes  is  in  order.  The  remainder  of  this  thesis  will  investigate  the  efficiency 
of  interleaved  RS  codes.  The  question  of  whether  or  not  the  efficiency  of  a  RS  coding 
scheme  increases  with  the  use  of  interleaving  will  be  answered. 

7. 1  Block  Interleaved  Guaranteed  Burst  Error  Correction 

Recall  that  block  interleaving  scrambles  codewords  to  be  transmitted  over  the  com¬ 
munications  channel  thus  randomizing  the  effects  of  a  burst  error.  Essentially,  the  burst 
error  is  broken  into  smaller  and  more  manageable  errors  over  multiple  codewords.  It  is 
the  depth  (d)  of  the  block  interleaver  that  determines  how  burst  errors  will  be  randomized 
over  the  interleaved  codewords.  Since  the  depth  of  the  block  interleaver  is  equal  to  the 
number  of  codewords  being  interleaved,  a  burst  of  length  b  =  rdm  bits  (where  r  =  2,3,4,... 
is  the  number  of  columns  in  the  interleaver  with  errors,  see  Figures  5.2  and  5.3)  will  pro¬ 
duce  burst  errors  of  b  =  rm  bits  maximum  in  each  de-interleaved  codeword.  Therefore, 
codewords  with  the  ability  to  correct  burst  errors  of  length  b  =  rm  bits  and  an  interleaver 
of  depth  d  is  all  that  is  required  to  correct  burst  of  length  b  =  rdm  bits. 

Reed-Solomon  codes  that  can  correct  t  =  r  random  symbols  errors  can  be  used  to  cor¬ 
rect  burst  errors  up  to  b  =  rm  bits.  However,  there  are  instances  when  burst  errors  do  not 
occur  exactly  over  t  =  r  symbols  and  the  random  error-correcting  capability  t  is  exceeded 
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(t  <  r).  The  equations  for  guaranteed  burst  error  correction  were  given  in  Section  4  and 
are  restated  here.  Guaranteed  single  burst  error  correction  is  b,  =  m(t  - 1)  +  1  bits,  guar¬ 
anteed  double  burst  error  correction  is  bj  =  m(Lt/2j  - 1)  +  1  bits  and  guaranteed  triple 
burst  error  correction  is  bj  =  ni(Lt/3j  - 1)  +  1  bits.  For  the  general  case  of  s  burst  errors 
the  equation  is  b^  =  ni(Lt/sJ  - 1)  +  1 .  These  equations  again  give  the  burst  error  correction 
of  a  given  RS  code  no  matter  where  the  burst  error  occurs  in  the  codeword. 

Combining  the  operation  of  the  block  interleaver  with  the  guaranteed  burst  error- 
correcting  capability  of  a  given  RS  code,  the  following  equations  are  derived. 

=  m{td  -  '\)  +  \ 


These  equations  express  the  single,  double,  and  triple  guaranteed  burst  error-correcting 
capability  of  a  block  interleaved  Reed-Solomon  code,  where  \\i^  is  the  general  form  of  the 
equation.  These  equations  simply  mean  that  using  a  knovra  RS  code  (n,k,t,m)  and  a 
block  interleaver  of  knoAvn  depth  (d),  a  burst  error  of  length  \\i  is  guaranteed  to  be  cor¬ 
rected  no  matter  where  it  occurs  over  the  interleaved  codewords. 
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12  Block  Interleaved  Guaranteed  Burst  Error  Correction  Efficiency 

In  Section  6  it  was  stated  that  from  the  Reiger  bound  and  the  Gallager  bound  ,  the 
phased  burst  error  correction  efficiency  of  a  linear  block  code  is  defined  as  z  =  2b/(n-k). 

If  z  =  1,  then  the  linear  block  code  is  said  to  be  optimal  giving  the  most  burst  error  cor¬ 
rection  for  the  amount  of  parity.  This  definition  of  efficiency  was  for  binary  linear  block 
codes  which  was  expanded  for  non-binary  RS  codes  later  in  Section  6  of  this  thesis 
(see  y, ,  72,  73 ,  and  ,  in  Section  6). 

Taking  the  definition  of  efficiency  (z)  another  step  forward,  it  will  now  be  expanded  to 
give  the  efficiency  of  block  interleaved  RS  codes.  Since  the  block  interleaver  increases 
the  number  of  parity  bits  being  transmitted  to  achieve  the  guaranteed  burst  error  correc¬ 
tion  given  in  \|/,  the  depth  (d)  is  applied  to  the  denominator  of  z.  As  before  m  is  applied 
to  the  denominator  because  RS  codes  are  non-binary.  Replacing  the  burst  error  correc¬ 
tion  factor  b  in  the  numerator  of  z  with  \\),  the  equations  for  the  block  interleaved  guaran¬ 
teed  burst  error  correction  efficiency  of  a  RS  code  have  been  derived  and  are  listed  below. 

T  = 

'  in(n  -  klyi 

4 

^  m{n  -  k)d 


m(n  -  k)d 
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The  importance  of  the  above  efficiency  equations  can  be  summarized  best  by  exam¬ 
ining  Figure  7.0  and  Table  7.0  below.  Figure  7.0  is  a  plot  of  the  block  interleaved  guar¬ 
anteed  single  burst-error  correction  efficiency  (  vg, )  of  length  255  for  RS  codes  of  vary¬ 


ing  interleaving  depths.  Length  255  is  chosen  here  because  of  the  popular  use  of  the 
eight  bit  (m  =  8 )  symbols.  Also,  only  the  high  code  rate  (R  >  75%)  RS  codes  are  shown 
to  reduce  the  congestion  of  the  graph. 

It  is  immediately  obvious  from  Figure  7.0  that  an  increase  in  interleaver  depth  pro¬ 
duces  an  increase  in  efficiency  of  any  given  RS  code.  Even  more  importantly,  examining 
the  data  from  Table  7.0,  it  is  seen  that  a  lower  code  rate  ( R )  non-interleaved  code  can  be 
replaced  (based  on  efficiency)  with  a  higher  code  rate  ( R  )  interleaved  RS  code.  This 
means  that  a  simpler  RS  code  can  be  interleaved  and  produce  the  same  burst  error  correc¬ 
tion  and  burst  error  correction  efficiency  as  that  of  a  more  powerful  non-interleaved  RS 
code.  By  using  a  simpler  RS  code  there  is  a  reduction  in  hardware  cost  and  the  number 
of  computations  needed  for  encoding  and  decoding  is  reduced.  As  a  reminder,  interleav¬ 
ing  will  produce  delays  in  communication  and/or  data  recording  systems.  Therefore,  re¬ 
placing  a  more  complex  RS  code  with  a  simpler  interleaved  code  can  only  be  done  if  the 
system  can  tolerate  the  interleaving  delays  introduced. 
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Table  7.0  Block  Interleaved  Guaranteed  Single  Burst  Error  Correction 
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BLOCK  INTERLEAVED  REED-SOLOMON  CODE 
Guaranteed  Single  Burst  Error  Correction  Efficiency 
(Code  Length  n  =  255) 


Figure  7.0  Block  Interleaved  Guaranteed  Single  Burst  Error  Correction  Efficiency 


7.3  Convolutional  Interleaved  Guaranteed  Burst  Error  Correction 

Convolutional  interleaving  is  very  much  like  block  interleaving  in  that  the  idea  is  to 
scramble  codewords  and  make  them  appear  random  to  the  communications  channel. 
Convolutional  interleaving  was  shown  in  Section  5  to  be  more  efficient  than  block  inter¬ 
leaving,  in  that  it  requires  less  storage.  For  the  convolutional  interleaver  it  is  the  depth 

(d)  and  the  number  of  storage  celts  (  ^  )  on  each  level  that  determine  how  burst  errors  will 
be  randomized.  From  Figure  5.4  it  can  be  seen  that  each  symbol  of  a  burst  error  wall  be 
separated  by  d  ^  symbols  after  convolutional  interleaving. 
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If  X  =  (Xo,  X,,  X2,  X3,  X4, ...)  is  a  codeword  applied  to  the  convolutional  interleaver  of 
Figure  5.4,  then  the  data  stream  to  be  transmitted  over  the  communications  channel 
would  be  that  of  Figure  7. 1 .  From  Figure  7. 1  it  is  clear  that  a  single  burst  of  length 


dC  symbols 

1 

d£  symbols 

n 

symbols 

EZ] 

dC  symbols 

ED 

m(df  +  1)  bits 


Figure  7.1  Output  of  a  Convolutional  Interleaver 


b  =  m(df  +  1)  bits  will  produce  a  single  random  error  in  each  RS  codeword  after  deinter¬ 
leaving.  Combining  this  fact  about  convolutional  interleaving  with  the  concept  of  guar¬ 
anteed  burst  error  correction  (see  guaranteed  burst  error  correction  capability  in  Section 
6)  the  following  equations  are  derived.  These  equations  express  the  single,  double,  and 
triple  guaranteed  burst  error  correction  capability  of  a  convolutional  interleaved  RS  code 
no  matter  where  the  burst  error  occurs,  where  (|)s  is  the  general  form  of  the  equation  for  s 
burst  errors. 


=  n^{dl  + 1)  - 1]  + 1 


=m 


+  l)-l 


=m 


(df  +  l)-l 


+  1 


+  1 


+  1 
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7.4  Convolutional  Interleaved  Guaranteed  Burst  Error  Correction  Efficiency 


As  stated  in  Section  6  the  phased  burst  error  correction  efficiency  is  classically  defined 
as  z  =  2b/(n-k).  Throughout  this  thesis  the  classical  definition  of  efficiency  has  been 
modified  for  various  applications  of  RS  codes.  In  Section  6  the  efficiency  (z)  was  ex¬ 
panded  for  the  case  of  guaranteed  burst  error  correction  for  RS  codes  (y,,  y2,  yj,)- 

this  section  the  definition  of  efficiency  (z)  has  been  expanded  to  include  the  case  of  block 
interleaved  guaranteed  burst  error  correction  efficiency  of  RS  codes.  Similarly,  in  this 
section  the  definition  of  efficiency  (z)  is  now  expanded  to  include  the  case  of  convolu¬ 
tional  interleaved  guaranteed  burst  error  correction  efficiency.  The  equations  provided 
below  are  derived  by  substituting  (j),,  <t)2,  <|)3,  and  (j)^  for  b  in  the  equation  z  =  2b/(n-k),  m  is 
applied  to  the  denominator  because  RS  codes  are  non-binary.  The  term  d  +  1  is  applied 
to  the  denominator  because  there  must  be  d  ^  +  1  codewords  processed  before  achieving 
the  guaranteed  burst  error  correction  of  <|),,  <j)2,  (1)3,  and  (j)s .  These  equations  represent  the 
guaranteed  burst  error  correction  efficiency  achievable  using  convolutional  interleaving. 
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«>,  = 


2m^{dt  +  l)-  l]+  1 
m{n  —  k^dl  +  l) 


Am 


03  = 
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1)-1 


+  1 


m{n  -  k^idt  +  l) 


6m 


4>,  = 


{dl  +  \)- 


m(n  -  k^di  +  l) 


2sm 


O,  = 


+  1)-1 


+  1 


m(rt  -  k)(di  +  l) 


As  can  be  observed  in  Figures  7.1  through  7.4  and  Tables  7.1  through  7.4  the  effi¬ 
ciency  of  any  given  RS  code  increases  with  an  increase  in  convolutional  interleaving 
depth.  Also,  it  is  clear  that  as  the  additional  delay  per  each  set  of  registers  (£)  is  increased 
the  single  burst  error-correcting  efficiency  is  increased.  Most  importantly,  notice  that  as 
the  amount  of  convolutional  interleaving  (delay)  is  increased  (here  as  £  increases  in  fig¬ 
ures  7.1  to  7.4)  the  efficiency  is  very  high  for  the  simplest  RS  codes.  For  example  the 
(255,253)  RS  code  with  convolutional  interleaving  of  d  =  2  and  f  =  4  (figure  7.4)  has  a 
guaranteed  single  burst  error  correction  efficiency  of  more  than  ninety  percent  and  can 
correct  a  single  burst  of  sixty-five  bits  with  a  code  rate  (R)  greater  than  ninety-nine  per¬ 
cent.  Very  impressive! 
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Table  7.1  Convolutional  Interleaved  Guaranteed  Single  Burst  Error  Correction 
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Burst  Error  Correction  Efficiency 


(X)N\/C)LLmONAL  IMTERLBWED  RED-SCDLOMON  CXX3E 
Guaranteed  Single  Burst  Emor  Connection  Efficiency 
(Code  Length  n= 255,  1  =  1) 


(255,253)  (255,247)  (255,241)  (255,235)  (255,229)  (255,223)  (255,217)  (255,211)  (255,205)  (255,199) 


Reed^omon  Code  (n,!^ 

Figure  7.1  Convolutional  Interleaved  Guaranteed  Single  Burst  Error  Correction 

Efficiency  (^=1) 
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Table  7.2  Convolutional  Interleaved  Guaranteed  Single  Burst  Error  Correction 

Eff.  (A.=2) 
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CONVOLUTIONAL  INTERLEAVED  REED-SOLOMON  CODE 
Guaranteed  Single  Burst  Error  Correction  Efficiency 
(Code  Length  n  =  255,  1  =  2) 
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Figure  7.2  Convolutional  Interleaved  Guaranteed  Single  Burst  Error  Correction 

Efficiency  (£=2) 
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Table  7.3  Convolutional  Interleaved  Guaranteed  Single  Burst  Error  Correction 

Eff.  (£=3) 
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CONVOLUTIONAL  INTERLEAVED  REED-SOLOMON  CODE 
Guaranteed  Single  Burst  Error  Correction  Efficiency 
(Code  Length  n  =  255,  1=3) 


Reed'Solomon  Code  (n,k) 


Figure  7.3  Convolutional  Interleaved  Guaranteed  Single  Burst  Error  Correction 

Efficiency  (^=3) 
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Table  7.4  Convolutional  Interleaved  Guaranteed  Single  Burst  Error  Correction 

Eff.  (/=4) 
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(255,253)  (255,247)  (255,241)  (255,235)  (255,229)  (255,223)  (255,217)  (255,211)  (255,205)  (255,199) 

Reed-Solomon  Code  (n,k) 


Figure  7.4  Convolutional  Interleaved  Guaranteed  Single  Burst  Error  Correction 

Efficiency  (f=4) 


8  Conclusion 


Burst  errors  are  commonplace  in  today’s  digital  communication  and  recording  sys¬ 
tems.  Reed-Solomon  codes  are  known  to  be  the  best  at  combating  burst  errors  because 
they  are  non-binary  Maximum  Distance  Separable  (MDS)  codes.  It  was  the  intent  of  this 
research  effort  to  explore  the  burst  error-correcting  capabilities  and  efficiencies  of  Reed- 
Solomon  codes. 

This  thesis  discussed  the  phased  burst  error-correcting  capabilities  of  RS  codes.  The 
concept  of  guaranteed  burst  error  correction  was  then  examined.  The  idea  of  guaranteed 
burst  error  correction  was  then  expanded  to  include  the  popular  block  and  convolutional 
interleaving  techniques  illustrating  the  increase  in  guaranteed  burst  error  correction 
achievable  using  interleaving.  The  classical  definition  of  error  correction  efficiency  was 
then  expanded  to  include  the  guaranteed  burst  error  correction  of  RS  codes  and  inter¬ 
leaved  RS  codes. 

It  was  shown  that  RS  codes  are  very  efficient  at  correcting  burst  errors.  In  fact  many 
RS  codes  approach  optimal  efficiencies  while  correcting  considerably  large  burst  errors. 
All  of  the  equations,  tables,  and  graphs  generated  in  this  thesis  can  be  used  by  code  de¬ 
signers  as  additional  criteria  in  choosing  the  most  efficient  burst  error-correcting  system 
for  their  needs. 
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9  Appendix  A 


0.714285714 


BURST  ERROR  CORRECTION  (bits) 


n,k 


63,61 


63,59 


63,55 


63,53 
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6.  Conclusions  And  Plans  For  Future  Activity 

Program  participants  have  remarked  frequently  on  the  success  of  the  RAPCEval  contract  and 
on  the  quality  of  the  results  that  have  been  produced  to  date.  The  eooperation  among  various 
engineering  communities  and  their  fruitful  interaction  has  been  a  very  satisfying  byproduct  of  the  entire 
project.  Students,  university  feculty,  government  and  private  sector  engineers  have  united  in  the 
common  attack  on  a  number  of  priority  Air  Force  electronie  warfere  concerns.  The  student  master  s 
degree  research  has  fiicused  on  topics  that  have  immediate  mpact  on  the  transitioning  of  new  and 
inproved  software  and  hardware  technologies  into  fielded  systems. 

Recently,  several  students  have  had  the  privilege  of  presenting  their  work  at  national  and 
international  conferenees,  where  comments  have  been  quite  fevorable.  Ron  Brinkley  presented  a  paper 
entitled  “Reed-Soloman  Coding  and  OMSK  Modulation  for  Cellular  Digital  Packet  Systems”  at  an 
TFF.F.  Pacific  Rim  Conference  on  Communications,  Conputers,  and  Signal  Processing  (PACRIM’  97). 
Henderson  Beigamin  gave  a  talk  entitled  “Seleetion  of  the  Most  ‘EflBeient’  Reed-Solomon  Code  for  a 
Specifie  Application  Using  Neural  Networks”  at  the  International  Conference  on  Communieations 
(ICC  ’99),  in  Vancouver,  British  Columbia.  These  aceonplishments  are  appreciated  all  the  more  when 
it  has  been  noted  that  it  is  very  unusual  for  a  master’s  level  graduate  student  to  present  a  paper  at 
national  and  international  conferences. 

Twelve  master's  degrees  have  been  awarded  at  the  time  of  this  report.  The  title  of  each  report 
is  listed  in  the  table  together  with  availability  information  for  each  document. 
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At  the  present  time,  several  students  are  planning  to  start  new  research.  Reports  fiom  various 
colleagues  in  the  military  and  industry  indicate  that  additional  students  will  be  interested  in  participating 
in  projects  of  the  kind  sponsored  by  the  RAPCEval  program. 
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